QUANTUM ERGODIC RESTRICTION THEOREMS, I: INTERIOR 
HYPERSURFACES IN DOMAINS WITH ERGODIC BILLIARDS 

JOHN A. TOTH AND STEVE ZELDITCH 
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5_^ ' Abstract. Quantum ergodic restriction (QER) is the problem of finding conditions on a 

pi |. hypersurface H so that restrictions ^j\h to H of A-eigenfunctions of Riemannian manifolds 

•^r ' (M, g) with ergodic geodesic flow are quantum ergodic on H. We prove two kinds of re- 

sults: First (i) for any smooth hypersurface H, the Cauchy data {(pj\H,dipj\H) is quantum 
ergodic if the Dirichlet and Neumann data are weighted appropriately. Secondly (ii) we give 
conditions on H so that the Dirichlet (or Neumann) data is individually quantum ergodic. 
The condition involves the almost nowhere equality of left and right Poincare maps for H. 
The proof involves two further novel results: (iii) a local Weyl law for boundary traces of 
eigenfunctions, and (iv) an 'almost-orthogonality' result for Fourier integral operators whose 
(~| , canonical relations almost nowhere commute with the geodesic flow. 
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L Introduction 

This is the first in a series of articles on "quantum ergodic restriction" (or QER) theorems. 
Quantum ergodicity (in the Riemannian setting) is the study of asymptotics of eigenfunctions 
of the Laplacian A^ on Riemannian manifolds (M, g) with ergodic geodesic flow. The QER 
question is whether restrictions of A^-eigenfunctions to hypersurfaces H G M are quantum 
ergodic along the hypersurface. More precisely, the aim is to find conditions on H so that 
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2 JOHN A. TOTH AND STEVE ZELDITCH 

the Dirichlet data ^x\h^ or Neumann data d^ipxln or full Cauchy data {ipx\H,dl^ipx\H) is 
quantum ergodic along H. Here, d^ is a fixed choice of unit normal to H. 

In this article, we consider the case of smooth hypersurfaces H G Q in piecewise smooth 
Euclidean domains Q C M^ with 'totally' ergodic billiards, i.e. such that all powers 
pk . ]^*gQ _i. B*dQ of the billiard map are ergodic. We first prove that in a suitably 
weighted sense, the Cauchy data is quantum ergodic for any hypersurface. This is a general 
phenomenon that holds for all Riemannian manifolds with boundary (further results are 
proved in |CHTZj ). We then consider the subtler question of when the Dirichlet data (or 
Neumann data) is itself quantum ergodic. It is simple to find examples of H for which (Pj\h 
are not quantum ergodic. But in Theorem H] we give a generic sufficient condition on H in 
terms of a 'dynamical' condition on the 'transfer map t^ : B*H — t- B*dVL from tangent 
CO- vectors to H to the boundary (see Definitions [2] and [3]). The proof is based on a seem- 
ingly new principle in quantum dynamics: the almost orthogonality, {F{Xj)ipj,ipj) — )■ 0, of 
A- eigenfunctions ipx (and their Cauchy data) to their images F{X)ipx under a semi-classical 
Fourier integral operator whose underlying symplectic map almost nowhere commutes with 
the geodesic fiow (or billiard map). 

To our knowledge, the only prior quantum ergodic restriction theorem concerns the Cauchy 
data along the boundary of eigenfunctions of the Laplacian on bounded domains satisfying 
Dirichlet or generalized Neumann boundary conditions |GL|,IHZ|,IBuj . However, heuristic and 
numerical work of Hejhal-Rackner jHR] appears to show that restrictions of eigenfunctions 
on finite area hyperbolic surfaces to closed horocycles are quantum ergodic. We restrict to 
Euclidean domains in this article because there are special techniques available for them and 
there are many interesting examples. At the end of the introduction, we discuss related work 
(and work in progress) on quantum ergodic restriction theorems in more general settings. 

To state our results, we introduce some notation (See ^HUl for a notational index). We 
consider the eigenvalue problem in a piecewise smooth bounded Euclidean domain Q C M"" 
with Neumann (resp. Dirichlet) boundary conditions: 

-Av^A, = A^v^A, in fi, 

(1) 
d,y(px, = (resp. ipx^ = 0) on dQ. 

Here, A is the Euclidean Laplacian and di, is the interior unit normal to Q. We denote by 
{(fx} an orthonormal basis of eigenfunctions of the boundary value problem corresponding 
to the eigenvalues Aq < A^ < X^- ■ ■ enumerated according to multiplicity. We denote by 
f3:B*dn^ B*dn the bilhard map of Q (see ^ }2lD . 

The quantum ergodic restriction problem for Dirichlet data is to determine the limits of 
the diagonal matrix elements 

pfiOp^ia)) := {Opx^{a)^xM^^>^M)LHH), Op^ia) € ^'JH)) (2) 

of (semi-classical) pseudo-differential operators Op^{a) along H. That is, a G 5*^/ iT*H x 
(0, Xq^]) lies in the class of zeroth-order semiclassical symbols with polyhomogeneous expan- 
sions in A~^ with A G [Aq, cxd) We recall that |DSj 

S°f{T*H X (0, X,']) := {a{s, r; A) G C^iT^H x (0, Aq ^]); a ~n^o E, ^A'^ 

(3) 

\dfdSis,r)\<C,,^Ar)-'-^'^}- 
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For simplicity, we denote the corresponding space of semiclassical pseudodifferential opera- 
tors by 

^l{H):=Op,{S'f{T*Hx[0,\,']). 

The sequence of Dirichlet data ipx^ln is called quantum ergodic on L'^{H) if almost all of 
the matrix elements ([2]) tend to the integral uj^{a) = J^,jjaod^ao of the principal symbol 
Qq against the quantum limit measure dfioo where d/ioo is the image under the transfer map 
T^ (see Definition |2]) to B*H of the limit measure d^B of |HZ] with boundary conditions 
B. The quantum ergodic restriction problem for Neumann data is similar, while that for 
Cauchy data is to study the matrix elements (A2CDh(v?Aj), CDh{'^Xj)), where A2 is a 2 x 2 
matrix of pseudo-differential operators on H. As can be seen in the table below, we need to 
weight the Dirichlet and Neumann data in a special way and to choose special A2 to obtain 
an almost unique weak limit. 

To put the quantum ergodic restriction problem into context, we recall that the Dirichlet 
or Neumann eigenf unctions are quantum ergodic on M (i.e. in the 'interior') whenever the 
geodesic flow of (M, g) (resp. billiard flow, if dM ^) is ergodic. In this article we restrict our 
attention to bounded domains in M", and refer to |GL[IZZ] for the relevant interior quantum 
ergodic theorems. It is shown in |GLt IHZt IBu] that (M, dM) satisfies QER if the classical 
billiard flow of M is ergodic, with u^ depending on the choice of boundary conditions (which 
we suppress in the notation). The boundary quantum ergodic limit measures are as follows 

m- 



Boundary traces and quantum limits 


B 


Bu 


u' 


d^B 


Dirichlet 


Uy 


X~^ dyU Y 


'f{q)da 


Neumann 


dyU Y 


Uy 


li^y'da 



Here, we denote by da = dqdrj the natural symplectic volume measure on B*dQ. The 
standard cotangent projection is denoted by vr : T*dQ — ?■ dQ. We also define the function 
'y{q,ri) on B*dil by 



'l-|r/|2. (4) 

The set S*M\h of unit vectors to M along H is a kind of cross-section to the geodesic 
or billiard flow and one might expect that on both the classical and quantum levels, the 
ambient ergodicity induces ergodicity on the cross section. However, this heuristic idea is 
based on the assumption that the Cauchy data (or Dirichlet data) on a hypersurface is a 
suitable notion of quantum cross section. In effect we are investigating the extent to which 
that is true (the abstract microlocal notion of quantum cross section is studied in |CHTZ2] ). 

1.1. Quantum ergodic restriction theorem for Cauchy data. The difference in the 
quantum ergodic restriction problems for Cauchy data versus Dirichlet data is visible in the 
one- dimensional case where M = S*^ or where M = [0, vr]. Modulo time- reversal, the classical 
systems are obviously ergodic and consist of just one orbit. The real- valued eigenfunctions 
sin kx, cos kx do not individually have quantum ergodic restrictions to a hypersurface (i.e. a 
point) X = xq. Indeed, they oscillate as fc — t- 00. On the other hand, if we weight the real 
Neumann data by ^ then the formula cos^ kx + sin^ kx = 1 shows that the Cauchy data is 
quantum ergodic. 
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A somewhat more revealing example is that of eigenfunctions of semi-classical Schrodinger 
operators h'^D'^ + V^ in one dimension and with connected level sets of the energy ^^ + V{x) = 
E. The real valued eigenfunctions have the WKB form Aj^ cos{ j^S{x)) and we see that again 
the Dirichlet data at a point fails to be quantum ergodic but that the Cauchy data is quantum 
ergodic if the Dirichlet data is weighted by S'{x) and the Neumann data by h. 

As these examples show, we should study the normalized Cauchy data 

CDH{ipx):MH,X-'d,ipx\H). (5) 

Henceforth, Cauchy data will always refer to this normalization. 

To state our result, we introduce some further notation. We assume that H G Q°, although 
with some extra technical complications, we could allow H to intersect dfl (e.g. the midline 
of the Bunimovich stadium). We denote hj s ^-^ qnis) & H a smooth parameterization 
of the hypersurface H G Q with s G f/ C ]R"~^ an open set and U E y ^-^ q{y) G dQ a 
local smooth parametrization of the boundary of Q. The dual coordinates are r G T*{H) 
and 1] G T*{dfl) respectively. We sometimes abuse notation and use (s, r) as coordinates on 
B*H, resp. {y,r]) as coordinates on B*dQ. 

The examples above indicate that we can only expect QER for rather special matrices 
A2 of pseudo-differential operators. We say that A2{h) is a normalized scalar matrix of 
semi-classical pseudo-differential operators if, 

• (i) A2i{h) = Ai2{h) = 0; 

• (ii) aA^As,T) = (1 - \r\'^)(rA22is,r) for all (s, r) G B*H. 

Theorem 1. Assume that Q C M" is a piecewise smooth Euclidean plane domain with 
totally ergodic billiards, and let H C Q° be any smooth hypersurface. Also, let A be a 
normalized scalar matrix of zeroth order pseudo- differential operators on H . Then, there 
exists a density-one subset S of N such that the (normalized) Cauchy data of Dirichlet or 
Neumann eigenfunctions satisfies 



lim {A2CDHi(p\XCDHi(px,) = / (o-^h + CrA22)c?Aioo, 

where dfi^o = ijdn)*^l^B is the push-forward by Tqq (DefinitionWi) to B*H of the limit measure 
dfiB in the table above. 

In other words, there exists a measure d/Xoo on B*H so that along a subsequence of 
eigenvalues of density one we have, 

{Opx{{l - \T\^)a{s, T))ipx\H, <^x\h) + \''^{Opx{a{s, T))dyjjipx\H, d^H^xU) 

(6) 
-^ /b,^2(1 - |rp)a(iyUoo- 

When the correspondence Tq^ has two branches, so that T^h T h;iVJT H-a where Tq^-';j = 
1,2, are single- valued maps (see section [3721) . the push-forward is given by 

y'^an)* '■— van >* + i'^an >*■ 

Similarly if there are more intersection points. 

Although we do not do so here. Theorem [1] can also be proved by using interior quantum 
ergodicity ^ZZj and a Rellich-type commutator argument. As a consequence, we show in 
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[CHTZ2J that this theorem can be naturally interpreted as a quantum ergodicity result for 
the quantum flux norm of the Cauchy data along the hypersurface, H <Z fl. Moreover, the 
latter argument is insensitive to whether or not dQ ^ 0. Thus, our result in |CHTZ2] extends 
Theorem [1] to interior hypersurfaces if C M of arbitrary compact ergodic manifolds (M, g) 
with or without boundary. 

1.2. Quantum ergodic restriction for Dirichlet data. The one-dimensional case raises 
doubts that there can exist a quantum ergodic restriction theorem for the Dirichlet data 
alone. However, it is not hard to show that quantum ergodic restriction does hold for 
random waves on general smooth hypersurfaces in dimensions n >2. In this section, we give 
a sufficient condition on H so that QER is valid. 

To motivate the condition, we consider another special case of the QER problem: i.e. the 
case of Neumann eigenfunctions of a Bunimovich stadium S" C M^, with the curve H given 
by the vertical midline of S. We note that S has a left-right symmetry ai (and an up-down 
symmetry). The midline H is pointwise fixed by ai- 'Half of the Neumann eigenfunctions 
are even with respect to ul and half are odd (in the sense of spectral density). Obviously, 
the odd eigenfunctions restrict to zero on H. On the other hand, the even ones are quantum 
ergodic on H: in the quotient domain S/Z2 by ai, H becomes a boundary component 
and even Neumann eigenfunctions of S are Neumann eigenfunctions of S/'L2- Hence a full 
density subsequence of their restrictions to H have the boundary quantum limit d^B given 
in the table above. Thus, there exist two subsequences of eigenfunctions of density 1/2 with 
different quantum limits, and therefore QER (quantum ergodic restriction) of the Dirichlet 
fails for (5*, H). Consistently with Theorem [H the Cauchy data is quantum ergodic. Indeed, 
the Neumann data is quantum ergodic on the midline by the arguments of |HZ[ IBuj applied 
to the half-stadium with Neumann boundary conditions on all but the midline and with 
Dirichlet boundary conditions on the midline. A more general example of the same kind is a 
curve H which intersects the midline in a segment of the midline. The same kind of example 
exists on the fixed point set of any manifold (with or without boundary) that possesses a Z2 
symmetry. 

We now give a sufficient condition for restricted quantum ergodicity of the Dirichlet data 
which, while difficult to verify in many cases, is generic for smooth hypersurfaces, at least 
for some classes of domains. As above, we let H G Q° (the interior of Q) denote a smooth 
hypersurface. The key objects associated to H are the transfer maps Tqq and r^. They 
are somewhat complicated because of the possible number of intersection points of a billiard 
trajectory with H. In §3.H we will give a more detailed definition of the transfer maps (see 
Definition 13. 5p . 

Let {y,ri) G B*dQ and let q{y) + t({y, rj) denote its billiard trajectory where (see §2.ip 

C{y,v) ■=V + V'^- W^y (7) 

is the lift of (y, 77) to an interior unit vector to Vt at y. 

Definition 2. The transfer map Tq^ : B*dVl — > B*H is defined as follows: We say 
that {y,ri) is in the domain of Tq^, (ie. {y,ri) G T^{Jqq)) if its billiard trajectory inter- 
sects H. Denote by E{tH{y,ri),y,r]) G H any intersection point (see section l37^) and by 
'^E(t (v»?) vr))C(l/)''?) th^ tangential projection of the terminal velocity at E(tH{y,ri),y,ri) to 
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TE{tjj{y,r]),y,ri)H . The graph of transfer map is by definition the canonical relation, 
Tr^ := {(2/,^,s(2/,r/),7r|(,^(y_^)_j^,,)C(2/,r/)) G B*dn x B*H, 

iy,v) e viri,), qH{s{y,v)) = EitHiy,v),y,v)}- (8) 

In the inverse direction, we have the double-valued transfer map (or correspondence) 

T^^ : B*H -^ B*dQ, (9) 

defined by taking {s,t) G B*H to the two unit covectors ^±{s,t) G S*,g.fi which project to 
it, following each of their trajectories until they hit the boundary and then projecting each 
terminal velocity vector to B*dQ. The graph of rfp is the canonical relation, 

On := {(s,r,y(s,r),7rj . ))C(s,r)) G B*H x B*dn}. (10) 



'H 



We emphasize that rfp is double-valued due to the fact that each (s, r) G B*H lifts to 
two unit covectors 

e±(s, r) = T+ Vl-|r|2zy± G Sl^ (11) 

which project to r. More precisely, the normal bundle N*H to an orientable hypersurface H 
decomposes into two R+ bundles N^, which we view as the two infinitesimal sides of N*H. 
Then ^±(s,r) G N^H are the lifts of r to unit covectors on the two sides. The lifts are 
analogous to (JTj) but only the interior side of the boundary is relevant in ([71) . The two lifts 
in (ITTl) give rise to two branches 



r± : B*H -^ B*dn (12) 

of r^^ by taking t±{s, t) to be the element of t^^(s, r) coming from the trajectory defined 
by the unit vector ^±{s,t) G S*^(^^^Q. This double- valued-ness is independent of the shape 
of H and is an essential feature of the quantum restriction problem. 

On the other hand, r^ is multi-valued due to the shape of H, specifically to the many 
possible intersection points of a billiard trajectory with H. This is an unavoidable feature 
of the QER problem which causes notational inconveniences more than essential analytical 
problems. For expository simplicity, we will assume that H is weakly convex, so that a 
generic line intersects H in at most two points. This saves us from tedious notations and 
indices for the branches of Tq^ caused by the multiple intersection points. In the course of 
the proof we indicate the modifications necessary to deal with general smooth hypersurfaces. 
In Definition 13. 5[ we introduce notation for the multiple branches. 

Another inevitable and tedious complication is that the correspondences have singular- 
ities at directions where billiard trajectories intersect H or dQ tangentially or which run 
into a corner of dQ. It is not apriori clear that such tangential rays cause only technical 
complications. Without a detailed analysis, it is possible that eigenfunction mass could get 
concentrated microlocally in tangential directions to H. In the case H = dQ it was proved in 
|HZ] that no such concentration occurs. The same method shows that no such concentration 
occurs on H either. In §8.41 we give a self-contained pointwise Weyl law argument to show 
that for a full density of u^'s, no such tangential mass concentration occurs. This method 
works equally well for manifolds, M, with or without boundary. 
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Another important symplectic correspondence associated to H is the once-broken trans- 
mission billiard map (or correspondence) through H, Ph '■ B*dQ — )■ B*dQ, defined by 



T-T+ {y,v), {y,r])e range r. 



+ 5 



l3H{y,v) = { (13) 

T+T-^iy,v), iy,v) e range r_. 

Equivalently, Ph rnaps the two endpoints of the rays through C,±{s,t) with (s,t) G B*H to 
each other. Thus, Ph follows the trajectory of {y, t]) G B*dQ backward to its intersection 
point qH{s{y,ri)) G H (assuming it intersects H), breaks at qH{s{y,ri)) by the law of equal 
angles, and then proceeds along the second link in the forward direction to dfl. Since 
the trajectory defined by (s, r) G B*H may have multiple intersection points with H, this 
map is also a multi-valued correspondence. With the simplifying assumption that H is 
weakly convex, a generic trajectory intersects H in at most two points, and we denote the 
corresponding branches of Ph by /3|^; k = 1,2. See §4.11 for further details. Note that the 
trajectory breaks only once even if the trajectory intersects H twice. 

Remark: Ph is not the same as the broken billiard fiow with a break on H because the 
trajectory travels backwards on the first link. Starting at {y,ri) G B*dQ, the map first 
time reverses the rj — )■ —77, then proceeds on the broken trajectory and then projects the 
terminal velocity. Hence, PH{y,v) differs by a time- reversal from the usual broken billiard 
trajectory. As a check, let us note that in the case of the midline of the stadium, the two 
unit vectors projecting to f G B*H are related by ai- Hence their rays, and the projections 
of the terminal velocity vectors, are related by ai- But the two-link ray starting at f G TydQ 
between them is not so related: the terminal velocity is the time reversal of crL{v). 

Given A C B*dQ in the following we denote the symplectic volume \dydr]\- measure of A 
by \A\. 

Definition 3. We say that {Ph,P) almost nowhere commute if the (commutation) sets 

CO,,, := {iy,v) e B*dn : PnP'iy^v) = P^Pniy^v)} 
have symplectic volume measure zero, \COp,k\ = 0, for all p,k = 1,2,3, .. . 

The almost nowhere commutativity (ANC) condition may be re-formulated in terms of 
'left' and 'right' return maps to B*H. Given covector (s, r) G B*H, one may lift it in two 
ways ( ITTj) to a unit covector (s,^±) G S}jQ which projects to r. We then follow the two 
geodesies 7^(t) with initial conditions (s, ^±) until they return to H on the same side as (s, ^±) 
and then project the terminal co-vector back to B*H. The almost nowhere commutativity 
condition is equivalent to the condition that these two oriented return maps do not coincide 
on a set of positive measure in B*dQ. The equivalence of the definitions is easily seen by 
re-writing the equation PHP''{y,v) = f^^f^H{y,v) (see ( 1T^ ) as 

T.T-'p\y,7^) = pPT^r-\y,7^) ^^ t-^P\+{s,t) = tZ^ PPt.{s,t), 

with (|/,?7) = Tj^{s,t). We observe that t^^P^t±{s,t) is the ± — > ± sided return map to 
B*H after k bounces off dVL. The union over k gives the full one-sided return map. The fact 
that k and p may differ in Definition [3] then seems rather natural, since the equality of the 
one-sided return maps should not depend on the number of bounces off dVt. 
Our main result on Dirichlet data is: 
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Theorem 4. Let Q C M" be a piecewise- smooth billiard with totally ergodic billiard flow 
and let H C int{Q) be a smooth interior hypersurface. Let ^\j]i = 1,2,... denote the L^- 
normalized Neumann eigenfunctions in fl. Then, if {f3H,(3) almost nowhere commute, there 
exists a density-one subset S ofN such that for Aq > and a{s, r; A) G Sj {T*H x (0, A(^"^]), 



lim {Opx^{ao)<^x,\H,Vx,\H)LHH} 



Aj— s>oo;jg5 

Here, Cn = yoi{S"-^) voi{n) 

of Dirichlet eigenfunctions, except that '-/{y, r]) 



c„ / a{s,T) pQ^{s,T)dsdT. 

' B*H 



and pg^dsdr := (tJ^)*(7 ^dydrj). Similarly for Dirichlet data on H 

, and for Neumann data if one 



m 



-1/2 



multiplies the above by (1 — |r| ). 

We summarize the conclusion in the following table. 



H traces and quantum limits 
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u^ 


dfioo 


Dirichlet 


Dirichlet 


U H 


{^dQ)*l{<l)dcr 


Dirichlet 


Neumann 


X'^dygU H 


(1- r-^)ir,^^)Mq)d^ 


Neumann 


Dirichlet 


U H 


{ra''n)Mq)-'da 


Neumann 


Neumann 


\~^dygU H 


(1- r-'){T,^^)M<l)-'da 



There is a stronger but more technical version of this result which uses a quantitative 
refinement of the almost nowhere commuting condition (see the definition of quantitative 
almost nowhere commutation (QANC) condition in Definition [T3] of subsection 18. 2. 3p . 

Theorem 5. Theorem^is valid as long as {Ph,P) quantitatively almost nowhere commute. 

To illustrate the result, we reconsider the midline H of the Bunimovich stadium 5*. In 
this case, for {y,ri) G B*, -.H, the two unit vectors ^±{y,rj) e S*fi projecting to rj are 0"^- 
related, i.e. images of each other under ai- Hence their trajectories, and the projections to 
TdQ of their terminal velocities, are ai related. It follows that (3H{y,v) = (^Liy.'u)- Hence 
P^ai = o'lP'' for all fc, and (S*, H) violates the condition of quantitative almost nowhere 
vanishing. Thus, our condition rules out this 'counterexample' to QER and the original 
ANC condition. At this time, we are not aware of any other kind of 'counter-example'. 
Indeed, the equality of the left and right sided return maps is a kind of left /right symmetry 
condition on H . It might be equivalent to the existence of an involutive isometry a : M ^ M 
fixing a positive measure of H, but at this time we do not know how to prove that. 

Unfortunately, it is often difficult to check whether a given hypersurface H in the billiard 
setting satisfies the almost nowhere commutativity condition. The next result shows that 
the condition holds for generic H in a. reasonably broad billiard setting. In ^ we prove: 

Proposition 6. Let Q be a two-dimensional dispersing billiard table (i.e. with hyperbolic 
billiards). Then for generic convex curves, H G Q, {Ph, P) almost never commute. 

This result is just an example of the genericity of the condition. We only use hyperbolicity 
of the billiards to prevent the existence of a positive measure of self-conjugate directions, 
and we restrict to dimension two to simplify one step in the argument. A much more general 
result on genericity of the condition should be possible, and probably could be proved by 
the methods of |PS1[ IPS2] . The non- commutativity condition is also similar to the condition 
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that Poincare maps of periodic reflecting rays of generic domains do not have roots of unity 
as eigenvalues; see |PSltlPS2| [S]. But that would take us too far afield in this article and we 
only give the simpler result to convince the reader that our non-commutativity condition is 
far from vacuous. In the boundaryless case, we believe that the condition can be proved to 
hold for closed geodesies of hyperbolic surfaces which are not fixed by isometric involutions 
|TZ2] and possibly for closed horocycles. 

We also remark that the 'total ergodicity assumption' is not overly resrictive: all of the 
known examples of ergodic billiards that we are aware of are totally ergodic. These include: 
the Bunimovitch stadium, the cardioid and generic polygonal billiards |CFSl[P] . The first two 
examples are known to satisfy stronger mixing conditions (the stadium is in fact Bernuollii 
|BU] ). It is apparently open as to whether generic polygons are mixing, but Troubetzkoy 
^] has recently proved that these billiards are also totally ergodic. 



1.3. Convexity of H and QER. As mentioned above, we assume H is convex for simplicity 
of exposition, but the assumption is not necessary for the quantum ergodic restriction to hold. 
Indeed, we may localize (and microlocalize) the problem by considering symbols supported on 
small regions of H . On a sufficiently small region of B*H, concave or convex hypersurfaces are 
essentially the same. Flat regions are different, but since tangential intersections only occur 
for a zero measure of rays, it suffices to prove that only a sparse (density zero) subsequence 
of eigenfunctions can concentrate microlocally on tangential rays. We prove such a result in 
Lemma 18.51 using a pointwise local Weyl law. In |HZ] the analogous proof is in Section 7 
(see Step 2 and Lemma 7.1), Section 9 (see Lemma 9.2) and in Appendix 12 (see especially 
Lemma 12. 5 in the Dirichlet case). 

The fact that QER is independent of the curvature of the hypersurface may seem surprising 
in view of the results of Burq-Gerard-Tzetkov [BGT] (see also [Rl IKTZt ITotlSo]) on L^ norms 
of restrictions of eigenfunctions. In these results, the L^ norm of the restricted eigenfunctions 
depends on whether the curve (or hypersurface) is geodesically convex or not. There is a 
dramatic difference, for example, between norms of restrictions to distance circles and to 
closed geodesies. 

However, these results do not make use of the global assumption of ergodicity of the 
geodesic flow. They are sharp in the case of the round sphere, but need not be sharp 
for negatively curved surfaces (for instance). The results of |HZt IBuj amply indicate that 
geodesic curvature of H is irrelevant to quantum ergodic restriction theorems: Indeed, when 
H = dfl, QER is valid for non-convex H and for H with flat sides (e.g. the Bunimovich 
stadium). 

Since the validity of the results for general smooth hypersurfaces is possibly surprising, 
and since we assume convexity at some points to simplify the notation, we signal in each 
section where and when we use the convexity assumption. 

1.4. Applications. Under the assumptions of Theorem |U the L^-restriction bounds along 
if C fi of a quantum ergodic sequence of eigenfuncitons are uniformly bounded above and 
below. Indeed, choosing a(s, r) = 1 in Theorem H] gives the following 

Corollary 7. Under the same assumptions as in Theorem^ 

, lim J|'^A,|//||i2(j:^) = c„ / p^^{s,T)dsdT. 
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Except for possible exceptional sparse subsequences of eigenfunctions, in the ergodic case, 
the asymptotics in Corollary [7] substantially improve on some recent bounds of U norms of 
restrictions of eigenfunctions to submanifolds in the articles [Rl IBGTt ITot T£\ [Soj IBRj ). 

In subsection I8.5[ we give an application of Theorem H] to the study of asymptotic nodal 
structure of eigenfunctions of ergodic planar billiards fi C M^. Consider either Dirichlet 
or Neumann eigenfunctions (p\ with —A^ip\ = X'^ip\. In |TZj we show that for general 
piecewise-analytic domains, the eigenfunction nodal set M^^ = {x E fl](p\{x) = 0} (where 
we omit the boundary itself) for both Neumann and Dirichlet boundary conditions satisfies 

caTd{Af^^ndQ) = 0{X). (14) 

Here, flT4|) requires no dynamical assumptions. In Theorem 6 of |TZj we proved a similar 
bound for intersections of nodal lines with interior analytic curves under an additional tech- 
nical assumption of 'goodness' of the curve C (see the proof of the corollary). The quantum 
ergodicity result in Theorem H] yields the following 

Corollary 8. Let Q c M."^ be a piecewise- smooth, totally ergodic planar billiard and C G Q 
a real- analytic interior curve. Assume {Pc,P) almost never commute. Then, for a full- 
density of (Dirichlet or Neumann) eigenfunctions, 

card{Af^^nC) = 0{X). 

1.5. Sketch of proof. We now outline the proofs of Theorems [1] and HI This outline is an 
integral part of the proofs of the Theorems; later sections fill in the considerable number of 
technical details, but rigidly adhere to the following outline. 

The first step is to use the 'quantization' of the transfer map rj^ : B*dQ — > B*H to 
transfer the Dirichlet data u^ := (P\\h of the eigenfunctions on H to the Dirichlet data 
u\ := (px\dn on dQ. As in |HZ[ITZj ). the semi-classical quantization of the transfer map is a 
layer potential integral operator. By Green's formula, we have 



«f(te) = Wc(A)«^](te):= / N^ni^,qH,q)u',{q)da{q), (15) 

J an 

where q G dfl, qu E H and where Nq^{X) is a boundary integral operator defined in Definition 
13.11 and f pOj) . There is a similar operator Nq^^ that transfers the Dirichlet data on the 
boundary to Neumann data u^^ '^ = d^ipxln on H. By a modification of a similar argument 
in |HZ] (see ©, we show that the operators N^^^iX), Ng^'^iX) : C^idn) -^ C°^{H) are 
semiclassical Fourier Integral operators quantizing rj^ : B*dQ — )■ B*H. 

We then consider matrix elements (|2]) of semi-classical pseudo-differential operators on H, 
i.e. operators with symbols a G S^l (T*H). From ( !T5|) it follows that for any a G S^l (T*H), 

{Op,ia)u^,u^)LHH) = {N^nW*0p,ia)N^^{X)uiu\)LH9n). (16) 

Similarly for Neumann data. The formula ( !T6|) is an advantage to working in the boundary 
setting; it allows us to use results of |HZj to avoid many technicalities. 

At first sight, it appears that this transfer operator fully reduces our problem to the 
known results in [HZ] on dfl. However, this is not the case. The key point is that 
Ng^{X)*Opx{a)NQ^{X) is not a pseudo-differential operator. Rather, in §11 we show that 

N^^{XrOp,{a)N^^{X) = Opxiir^^ya ■ p^J + ^^(A), (17) 



QUANTUM ERGODIC RESTRICTION 11 

where pg^ G C°°{B*H) is a smooth density and where F2{X) is a semi-classical Fourier 
integral operator which quantizes the transmission map (fT3!) . Thus, one of the main contri- 
butions of the present article is the analysis of the matrix elements of -F2(A). 

The proof of Theorem [1] is then rather simple: The Neumann data produces a second 
Fourier integral operator F2^(A). The key observation is that the composition of F2{X) on 
the left by an operator with symbol 1 — jr^p cancels F2{X) in the sum of the matrix elements 
with respect to both the Dirichlet and Neumann data. Thus, the Fourier integral operator 
terms cancel, leaving only a pseudo-differential matrix element whose limit was calculated 
in [HZ] . 

For the Dirichlet data alone, there is nothing to cancel F2{X). Therefore, an important 
step of the proof of Theorem H] is to determine the weak* limits of the matrix elements 
{F2{X)u\,u\). This is a special case of the more general problem of determining limits of 
matrix elements of Fourier integral operators relative to eigenfunctions, which is discussed in 
|Z2] and will be investigated more thoroughly in |TZ3] . A significant additional complication 
is that the matrix elements are relative to boundary restrictions of eigenfunctions rather than 
to the eigenfunctions themselves. 

In the case of pseudo-differential matrix elements, the first step in obtaining the limits of 
matrix elements relative to A- eigenfunctions is to show that they are invariant measures 
for the billiard flow. In the case of matrix elements of Fourier integral operators, it is not 
clear what kind of invariance properties the limit has, viewed as a functional of the symbol 
of F2{X), when the canonical relation of F fails to be invariant under the billiard flow. Under 
the condition of Definition ([3]), the canonical relation is, in fact, almost nowhere invariant 
under the billiard flow. We turn this to our advantage by using a rather surprising averag- 
ing argument to show that the limits must be zero under the almost nowhere commuting 
condition. 

The starting point of this argument is the classical boundary integral equation 



<{q)= / N'enW{q,q')ui{q')da{q') (18) 

J on 

for the Dirichlet data of ipx on dfl. The operator Nq^{X) is studied in detail in [HZ] as a 
semi-classical Fourier integral operator quantizing the billiard map /3 : B*dQ — )■ B*dQ. The 
classical formula, Nq^{X)u\ = u\ was the key input to the quantum ergodicity result of [HZj . 
In particular, it shows that when H = dfl, the weak* limits of ([2]) are invariant under /3. 
We apply this quantum invariance to the -F2(A) term. For all m, we have 

(Op,(aX,«f)^.(H) = {Nl^{XrN^nW*Op,{a)N§,{X)NlSiXrului),^.^eny (19) 

At this point, we cannot apply Egorov's theorem to the conjugation because of the almost 
nowhere commuting condition in Definition [3l Under that condition, the canonical relation 
underlying Ng^{X)"^NQ^{X)*Opx{a)NQ^{X)NQ^{X)"^ changes with each m, and it appears 
that this destroys any invariance property of the quantum limits. What we show is that 
it forces the matrix element {F2{X)u\,u\) — )■ on average as A — )■ oo. Thus, F2{Xj)u\. is 
'almost orthogonal' to u\, under the almost nowhere commuting condition. 

The proof is based on a special case of a novel boundary local Weyl law for semi-classical 
Fourier integral operators. 
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Theorem 9. Let Q C W^ be a Euclidean domain with piecewise- smooth boundary and u\. = 
V^Aj|(9Q (resp. u\, = Y'di,ipXj\dn) be the Neumann (resp. Dirichlet) eigenfunction boundary 
traces. Let F{X) : C°°{dQ) — ?■ C°°{dQ) be a zeroth-order X-FIO with canonical relation 
Tf^x) ■■= graph{KF) C B*{dQ) x B*{dVL) and let 

Sf(a) := U {{y.v) e B*dn-r{y.v) = ^Fiy,v)}- 

Then, under the assumption that \Lp(\)\ = 0, 

^ ' A,<A 

In |TZ3] . more general local Weyl laws for semi-classical Fourier integral operators are 
proved. The general boundary local Weyl law is that there exists a density pp € C°°{B*dVt) 
such that 

tLWX^ S (^(^^K'^) = / <F{X)){y,r^) .pFiy,r^)dydv. (20) 

Moreover, ppdydrj can be explicitly computed as a sum of certain measures associated with 
the various multi-link contributions to the fixed-point set, Tip(^x)- The proof is valid for 
arbitrary Riemannian manifolds with piecewise-smooth boundaries. We present only the 
narrow version above since it is the only statement we need for Theorem H] and since the 
proof of the general result is lengthy. Indeed, the proof of the special case in Theorem [9] is 
already the lengthiest part of the proof of Theorem H] due to the complicated nature of the 
wave group on a domain with boundary. 

Using Theorem [9l we reduce Theorem H] to the following almost orthogonality theorem, 
which has an independent interest. 

Theorem 10. With the same notation as in Theorem\^ assume that F{X) is a zeroth order 
X-FIO and assume that its underlying symplectic map kf almost nowhere commutes with f3 
in the sense of Definition\^ Then 

-i- Y. l(^(AK.><.)r = 0(1) as A ^ oo. 

^ ' A,<A 

In Theorem 18.11 we restate the result in the form we need for the proof of Theorem |H 
However, the proof is valid with no essential modification in the general case. An interesting 
point is that the result does not use ergodicity of the billiard map. It is a completely 
general almost orthogonality result for pairs of operators satisfying the almost nowhere 
commutativity condition. 

The proof of Theorem[TO]is based on a somewhat curious averaging argument. By quantum 
invariance, we may (and do) replace F{X) by the operator 

1 ^' 
<^(^))^^ -■= M E[^l^(^)Ti^(A)[iV|^(A)]-. (21) 

m=0 
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Then, by the Schwartz inequahty, 

An^A A j ^ A 

(22) 
where, 

Fm(A) := (F(A));,(F(A))m. (23) 

Under the almost nowhere commuting assumption, these operators -Fm(A) satisfy the as- 
sumptions of Theorem [91 The almost nowhere commuting condition (and its quantitative 
refinements) arises because the fixed point set of the canonical relation of ( 123|) is 

M 

Sfm(a)= U [j^(^B*dQ; /3^/3UO = P8n/3\0} (24) 

k=-Mpez 

and the quantitative almost-nowhere commutativity condition (see Definition [T3ll implies 
that 

limsup— |S^^,(A)| =0. (25) 

M^oo ^yj- 

This is sufficient to imply that the average over eigenvalues \j of ( 12^ is 0{-^) and hence 
the almost orthogonality statement in Theorem [TOl 

Quite a lot might be lost in the use of the Schwartz inequality in (l23l) . As remarked in §8.21 
the almost nowhere commutativity condition is the condition under which averages over the 
spectrum of the norms 1 1 {F2k{Xj', a, £))m u^W tend to zero. It is another curious feature that 
we obtain a 'subsequence of density one' along which the almost orthogonality holds. It is 
not clear whether (as for quantum ergodicity) this is an essential feature or a defect of the 
proof. 

1.6. Final introductory remarks. We close the introduction by posing the quantum 
ergodic restriction problem in a more general context, and by describing some continuing 
work on the general problem. 

The quantum ergodic restriction problem concerns the linear functionals ([2]) defined on the 
space ^°c(if ) of zeroth order semi-classical pseudo-differential operators on H. As we have 
seen, there may exist large subsequences of eigenfunctions which restrict to zero; on the other 
hand, the known restriction estimates (see [BGT[[KTZj ) do not give uniform bounds on the 
L?'{H) norms of the restrictions. Hence the family {p^} ^ of functionals corresponding to 
Dirichlet (or Neumann) data of A- eigenfunctions along H is not apriori a bounded family 
of functionals, hence not necessarily compact in the weak* topology. Nevertheless, we may 
pose the 

Problem 11. Let H G M be a hyper surf ace. Determine the set Qh of 'restricted quantum 
limits', i.e. weak* limit points of the sequence of the (un-normalized) semi-positive states 
{p^} corresponding to the Dirichlet (or Neumann) data of eigenfunctions on H . Similarly, 
determine the restricted quantum limits for the functionals defined (as in TheoremU^ by the 
full Cauchy data ^. 

The analogue of an interior quantum ergodic sequence of eigenfunctions is given by: 
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Definition 12. The pair (M, H) satisfies QER (quantum ergodic restriction) for Dirichlet 
data if Qh contains a limit state u^ so that a full density subsequence of pf — )■ u^ in the 
weak* sense. Similarly for Neumann and full Cauchy data. 

In [TZ2j . we prove a quantum ergodic restriction theorem for a general Riemannian man- 
ifold (M, g) without boundary. In this setting, there do not exist transfer maps and broken 
billiard maps, and we base the proof on the almost nowhere equality of the left and right 
return maps to H. We believe the condition is satisfied for closed geodesies of negatively 
curved surfaces which are not fixed by isometric involutions. 

In [CHTZ] we use somewhat different methods from this article on manifolds with or 
without boundary to prove QER for the Cauchy data when the billiard flow is ergodic. 
We use two methods: one method (close to that of this article) to study the (matrix) 
Calderon projector C(A) : ip\ — )■ {^\\H)dyip\\H) as a semi-classical Fourier integral operator 
quantizing a certain first return map to S*fjM (the unit cotangent bundle of M along H). 
Our second method is a modification of the approach of |But IGLj in directly relating interior 
and restricted quantum ergodicity. 

In further work in progress |CHTZ2j we are studying a more general microlocal setting in 
which H (more precisely, S^M) is replaced by a general hypersurface T C S*M which is 
(roughly speaking) a cross-section to the geodesic flow. Quantum ergodicity of Cauchy data 
should hold if the first return map for S*^M is ergodic. The relevant notion of Cauchy data 
is provided by a microlocal restriction operator to T that arises in Grushin reductions. It is 
then an interesting question to relate this abstract reduction and the conditions for QER to 
those obtained in this article and in |CHTZj for the cross section defined by S^Vt, the unit 
cotangent vectors with footpoints on H . 

Acknowledgements We thank H. Christianson and H. Hezari for comments and collab- 
oration on related work. We also thank A. Strohmaier and S. Nonnenmacher for helpful 
comments on the averaging argument for Fourier integral operators. Thanks also to D. 
Hejhal for comments on the relation of our results to those of |HR] . and to Peter Sarnak for 
comments on QUE and our Corollary [HI 

2. Classical billiards on a domain with corners 

This section is a compilation of rather standard facts and definitions regarding billiards on 
domains with corners. The discussion does not involve H or any objects specific to our QER 
problem. However, we go over the basic definitions in some detail for the sake of clarity and 
to introduce notation that will be used in the rest of the article. 

We assume throughout that Vt C M'" is a smooth domain with corners. We define a 
smooth domain with corners and with M boundary faces (hypersurfaces) to be a set of 
the form {x G M™' : Pj < 0, j = 1, . . . , M}, where the defining functions pj are smooth in a 
neighborhood of f2, and are independent in the following sense: at each p such that Pj{p) = 
for some finite subset j G Ip of indices, the differentials dpj are linearly independent for all 
j E Ip. A boundary hypersurface Hj is the intersection of Vt with one of the hypersurfaces 
{pj = 0}. The boundary faces of codimension k are the components of Pji = ■ ■ ■ = Pj^ = 
for some subset {j'l, . . . ,jk} of the indices; each is a manifold with corners. This definition 
is extrinsic in the sense of [M] and is not the most general one possible. 
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It is essential to allow Q to have corners since domains with ergodic billiard flow in R" 
have corners. The singularities play a very significant role in the dynamics. They occur at 
corners, at trajectories which become tangent to the boundary, and in our setting of domains 
with obstacles, at trajectories which run tangent to the obstacle. 

We denote the smooth part of dQ by {dQ)°. Here, and throughout this article, we denote 
by W° the interior of a set W and, when no confusion is possible, we also use it to denote 
the regular set of dQ. Thus, dfl = {dfl)° U S, where E = [j^^jiHi fl Hj) is the singular set. 
When dimfi = 2, the singular set is a finite set of points and the Hi are smooth curves. In 
higher dimensions, the Hi are smooth hypersurfaces; Hi n Hj is a stratified smooth space of 
co-dimension one, and in particular E is of measure zero. We denote by S^Q the set of unit 
vectors to Q based at points of S. We also define C°°{dQ) to be the restriction of C°°(R") 
to 911. 

We define the open unit ball bundle B*{dQ)° to be the projection to T*dQ of the inward 
pointing unit vectors to Q along (dQ)". We leave it undefined at the singular points. 

To simplify the exposition, we describe B*dQ as if it consisted of tangent vectors rather 
than cotangent vectors, using the Euclidean metric to identify the two. We denote the 
standard projection by vr : B*dQ — )■ dQ. 

2.1. Billiard map. At smooth points of dQ, the billiard map (3 is defined on B*dQ as 
follows: given {y,r]) G B*dQ, i.e. with |?7| < 1 we let (y, C) ^ S*Q be the unique inward- 
pointing unit covector at y which projects to {y,ri) under the map Tq^Q — )■ T*dfl (see (j7])). 
Then we follow the geodesic (straight line) 

{y,r]) e B*{dnr ^ E{t,y,r]) := q{y) + tay,v) G fi C M" (26) 

to the first place it intersects the boundary again; let y' G dQ denote this first intersection. 
Thus, we find the least t = t{y,ri) so that pj{E(t,y,ri)) = for some j and let y' := E{y,ri) 
solve the equation 

(l{y') = E{t{y,r]),y,ri). 

The notation E(t, y, rf) is intended to suggest 'exponential map', i.e. the position in Vt reached 
by the billiard trajectory at time t. 

Denoting the inward unit normal vector at y' by Vyi, we let C' = C + 2(C ■ ^y')^y' be the 
direction of the geodesic after elastic refiection at y\ and let r/' be the projection of rj' to 
T*,dQ. Then we define explicitly, 

(3{y,r]) = iy',r]'). (27) 

We call each (linear) segment of a billiard trajectory a 'link', so that (3 relates the initial 
and terminal directions of a link. We denote the first link of the billiard trajectory determined 
hy{y,r])eB*dnhY 

{y,v)- (28) 

We also denote by 

G" : S*n -^ S*n (29) 

the (partially defined) broken geodesic fiow in the full domain 5**^2. The trajectory G'^(x,^) 
consists of the straight line motion between impacts with the boundary, and by the usual 
refiection law at the boundary. Thus, /3 is the reduction to the boundary of C^. We refer to 
|PSt ISV] for background. 
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We now address the complications caused by tangent directions and singular points. The 
billiard map is not well-defined on initial directions which are tangent to dQ, i.e. {y, t]) e 
S*dQ with |?7| = 1. If the domain is convex at y, then the line q{y) + trj immediately exits 
the domain. On the other hand, if it is concave at y, then the trajectory remains in the 
interior of Q and defines a link. If the domain is a 'saddle surface' at y, then some tangential 
links stay in Q and others exit. We introduce some terminology to distinguish these cases. 

Definition 2.1. Let {y,ri) g S*d^ be a unit tangent vector. 

(1) We say that {y, t]) is an exit direction, if the oriented ray q{y) + trj exits Q for t > 0. 
We denote the set of exit directions in B*dQ by £{y) and define [i{y,rf) = {y^rf) for 
such exit directions. 

(2) We say that {y, rj) is an interior tangent direction, if q{y) + trj remains in Q for 
sufficiently small t > 0. We define (3{y, rj) = {y, rj') by following q{y) + trj for t > 
until its first point of intersection y' G dVt and let rj' be the tangential projection of rj 
aty'. 

It is common to puncture out all tangent directions from B*dQ, i.e. all r] G S*dQ with 
|?7| = 1 and only define billiard maps on the open unit ball. But it makes sense to retain 
unit tangent directions as long as the straight-line trajectory remains well defined in Q. Unit 
tangent vectors to dQ at convex points can never be terminal velocity vectors of incoming 
trajectories and are therefore irrelevant to the definition of the billiard fiow. However, when 
dimfi > 3, as in the case where the boundary is a monkey saddle, it is possible for an 
incoming trajectory to arrive at an exit direction and terminate there. 

We will also need to distinguish initial directions whose trajectories terminate in tangent 
vectors to dQ. 

Definition 2.2. Let dimf2 > 3. We denote by Q{y) C B*{dVtY the set of grazing directions 
rj such that (3{y,r]) G S*dQ, i.e. such that the terminal vector of the link is tangent to dQ. 

No such exit directions exist when dimfi = 2. We now consider the regularity properties 
of P and the impact of the corner set. 

Definition 2.3. We denote by B^dfl or equivalently by S^Q the set of unit vectors to Q at 
points of the singular set. We define 

n^ = B* {{dny\E-\E))\g, 

i.e. the set of smooth directions of links whose next collision with dQ is neither at S nor 
tangent to dfl. 

Lemma 2.4. The domain of P is TZ^ U Q. It is a smooth symplectic map on the interior of 
TZ^ , which is smooth domain with corners in B*{dfl). 

Proof. It is clear from the definition that f3 is well-defined on 7?.^ U Q and on no larger 
domain. It is standard that /3 is symplectic where it is defined; see |CFSj . A simple way 
to prove it is that /3 has a generating function, namely the boundary distance function 
da{y,y') = \q{y) — q{y')\ on dQ x dQ. We have punctured out Q since /3 has jumps along 
this set, i.e. is discontinouous. 

Regarding regularity, it is clear from the formulae fl26l) - fl27j) that E{t,y,ri) and P{y,ri) are 
smooth wherever Uy is smooth and t{y, rj) is smooth. Of course, Uy is smooth away from the 
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corners. Clearly, t{y,ri) is smooth as long as Uy is smooth, E{t{y,ri),y,'r]) ^ E and as long 
as {y,ri) ^ Q; at the latter points, t{y,ri) is dis-continuous. Hence, /3{y,ri) is smooth on 71^ 
and on no larger set. 

We observe that E^^(T,) C B*dQ is locally defined as the simultaneous zero set of two 
smooth functions. Indeed, it the union of the sets {{y,ri) : E{y,ri) & Hi H Hj}. We may 
extend the hypersurfaces Hj beyond Q as the zero sets of the pj and consider the travel time 
tj from {y,ri) to these hypersurfaces. I.e. we may extend the time functions tj{y,7]) as the 
smallest value of t so that pj{E(t, y, r])) = 0. Then E~^{HinHj) is defined by ti{y, ri)tj{y, rj). 
Similarly for multiple intersections. 

D 

2.2. Iterates of the billiard map. We now consider the inverse /3~^ of (3 and their iterates. 
We note that the billiard fiow is 'time-reversal invariant', i.e. P~^{y,r]) = —(3{y, —t])- Hence 
/3~^ is well defined on the set 7^~^ of directions in B*{dQ)° such that E{y, —rf) ^ S. With a 
slight abuse of notation, we let /3~^(-B|;) = £'~^(S), i.e. it is the set of directions which map 
to the corner set; strictly speaking, there we do not define the image directions there. 

We now define open sets TV' such that TV D TZ'^ D . . . , TZ~^ D TZ~'^ D . . . where /3*^^) is 
defined, for any integer k. 

Definition 2.5. For k >1, we define the set TV' of k-regular directions such that the first 
k links of its trajectory stays in B*{dQ)°\Q . For k < —1, we define it similarly for (5~^ . 
The coball singular set is the set S* := [Jf^P~''{B^dQ) and the N-th order coball singular 

Definition 2.6. We define the generalized grazing set to be 

g* = [jp''{S*dQ). 



Similarily, the N-th grazing set is denoted by 

G*N= [j P\S*dn). 

\k\<N 

By repeating the proof of Lemma 12. 4[ we have: 
Lemma 2.7. The domain of [i^ is TZ^ It is a smooth symplectic map on the interior of TV. 

The regularity of 13^ near (y, r]) is essentially the issue of the regularity of the time t^^^ (y, r]) 
at which the billiard trajectory defined by (y, r/) hits dVt for the kth time. It fails to be defined 
if the trajectory runs into S or becomes an exit direction before this time, and additionally 
fails to be smooth if the trajectory becomes tangent to dVt at time t^''\y,r]). 

Since the corner set and its images and pre-images under /3 is a countable union of stratified 
smooth hypersurfaces, /3 is a measure-preserving map on a full-measure subset of B*{dQ)''. 

3. Hypersurfaces H: classical and quantum transfer maps 

In this section, we define the quantum transfer operators Ng^{X) (resp. Nq^) (see f lTSj) ). 
and show that they are semi-classical Fourier integral operators associated to the billiard 
relation given by transfer maps rj^ (resp. t^^). As discussed in the introduction, the 



Go{x,x',X) = ^A"-2(27rA|x-x'|)-("-2)/2Ha„/2-i(A|x-x'|) (31) 
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quantum transfer operators are fundamental in our approach, which is based on transferring 
the expected value {Op\{a)u^ , u^) l2(^h) to the boundary by the identity (TT6|) . 

The operator Ng^{X) : C°°{dQ) — )■ C°°{H) is defined as follows: By Green's theorem, if 
ip\. is a Laplace eigenfunction, then for any x G f2°, 

<^A,(a;)= / [d^yGo{x,q,\j)u\iq)-Go{x,q,\j)d^^ipxjiq)]da{q), (30) 

Jan 

where 

i 

4' 

is the free outgoing Green function on M". The Neumann (resp. Dirichlet) boundary con- 
ditions eliminate the second (resp. first) term, and by restricting z to the hypersurface in 
( 130|) . one gets the boundary integral equation 

The transfer operator from boundary data to Dirichlet data along H is defined as follows: 

Definition 3.1. Nq^{\) : C'^{d^) -)■ C°°{H) is the operator whose Schwartz kernel with 
respect to the hypersurface volume element da{q) on {OVty is given by, 

d„^Go{qH,q, Xj), Neumann 
N^ni<lH,q,X,) = { (32) 

Go{qH,q, Xj), Dirichlet 



Similarly, the transfer operator from boundary data to Neumann data along H is given 



m: 



Definition 3.2. With the same notation and assumptions, Ng^{\) : G°^{dn) -)■ G°°{H) 
is the operator with Schwartz kernel 

{d,.,^d^^Go{qH,q,Xj), Neumann 
(33) 
9u,^GoiqH,q,Xj), Dirichlet 

We now prove that Nq^{\j),Nq^{\j) are semi-classical Fourier integral operators (with 
singularities) quantizing the partially defined multi- valued transfer map Tqq between B*dVL 
and B*H. This is the first place where the convexity of H is relevant. 

We recall that a semiclassical (ie. A)-Fourier integral operator, -F(A) : G'^{dVL) — )■ G°^{dVL), 
is an oscillatory integral operator whose Schwartz kernel can be written locally in the form 

F(y, y'- A) = (27rA)--i /" e^^*(''(^)'^(^')'^)A(g(y), q{y'), 6, A) dO (34) 



n-l 



where the amplitude is a semi-classical (polyhomogeneous) symbol in A, A(g(y), g(|/'), 0] A) ~ 
Er=o M<li.y). Qiy'), ^)A"^-^ with Ak e G^{dQ xdQx M""^), A; = 0, 1, 2, .... The associated 
canonical relation is defined by 

Tp = {{y, Vy^, y, -Vy>^ : Ve$ = 0} C T*dn X T*dn. (35) 
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Proposition 3.3. Nq^{\j), resp. ^J^Ng^{\j), are semi-classical Fourier integral operators 
of order zero with canonical relation ^. The principal symbol of Nq^ is the half density on 
the graph of rf/^ given at a point {{y , v) , Tmiv , v)) ^ ^r^ with {y,T]) e B*dn by 



'an 



in the symplectic coordinates {y,T]) of B*dQ. Also, 



.(iVa''o(A))(y.,)=(-^^^)"V.*,| 



a(\-^N«„'(X)){y,ri) = i{^{y,rMrm{v,V))y'')dydri\''\ 

Remark: We observe that the symbol of Nqq{X) is infinite on the set where 7(tJ^(|/, tj)) = 0. 
This is also true in the case H = dfl in jHZj . We will cut off away from this set in the proof 
of the main theorems (see Proposition I4.2p . The key point is that for a quantum ergodic 
sequence, eigenfunction mass cannot concentrate in this set (see Lemma 18.51 in subsection 

I82D. 

Proof 

The proof of the Proposition consists of a sequence of Lemmas. We first observe that 
the quantum transfer maps are Fourier integral operators with phase equal to the distance 
function \qH — (l\ with qn & H,q E dQ. We then discuss the corresponding canonical relation. 
In section 13.31 (see Lemma 13. 6p we calculate the principal symbols. 

Under our assumption that H C Q°, we have \qH — q\ > dist{H, dQ) > 0, and we may use 
the asymptotics of the Hankel function flHT]) to obtain, 

oo 

Go{qH{s),q{yy,Xj) ~A,^oo (27rA,)^e^^^l^«(^)-^(^)l5^6fc(g^(s),g(i/))AT'= (36) 

fe=0 

with bk G C°°{H X dfl) and 6o > 0. The expansions for duyGo{qH{s),q{y); \j) and of 
9uy^di,yGo{qH{s),q{y); \j) are of the same form and are computed by taking the normal 
derivatives of the expansion in ( l36l) . Thus, the operators in Definitions (I3.ip -( l3l2l) are semi- 
classical Fourier integral operators with phase equal to the distance function \qH{s) — q{y)\ 
between a point of H and a point of dQ. The canonical relation is given by 

WF'.iN^^) C {(., <,(,)r(gH(3), qiy)), y, n^yfiqHis), qiy)))} C B*H X B*dn, (37) 

where, 

r{qH,q) ■= -. ^• 

m - q\ 

The canonical relation in ( 137|) is easily seen to be the same as (IHl). 

Remark: This statement is simpler than the one for the boundary integral operator Nq^{\) 
f lTSj) in [HZj. In that case, the phase and amplitude are singular on the diagonal and have to 
be regularized. Here, the operator kernel is smooth away from the corners of dQ. The singu- 
larities at the corners are only a technical inconvenience in the quantum ergodic theorems, 
and we refer to |HZ] for a more detailed discussion of them. 
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3.1. Details on the transfer map rj^. We now describe in more detail the canonical 
relation (137|) and the transfer map rj^ of Definition [2l As mentioned in the introduction, 
the definition of this map involves several technicalities: 

• The transfer map rj^ is only partially defined on B*dQ, namely on the range of the 
adjoint transfer map t^^ : B*H — )■ B*dVt. 

• Tqq is multi- valued since the trajectory {y, rj) defined by (y, rj) G B*dQ may intersect 
H generically in multiple points. 

• It has singularities when the trajectory {y,r]) from {y,r]) e B*dQ is tangent to H or 
when the trajectory is tangent to dQ. Such a trajectory may go on to intersect H 
again. 

• The adjoint transfer map r^^ : B*H — )■ B*dQ is singular when a ray from B*H runs 
into a corner. 

We now introduce notation to deal with these issues. At this point, we assume that H 
is convex. With this assumption, a generic line intersects H at most twice (indeed, it only 
intersects it more than twice if in intersects it in a fiat part). 

Definition 3.4. The domain C of Tq^ is the set of [y.-r]) E B*{dD.y such that the link 
{y, —rj) (i.e. E(t,y, —rj)) intersects H (see ( f^) ). Let C° d C denote the open dense subset 
such that (y, —rj) fl H consists of two distinct points {q\j{y, — r/), (fnijj, ""'?)}; ordered so that 

ki?(?/, ~v) -y\< \ql{y, -v) - y\, 

i.e. so that the line intersects H first at q}j{y, —rj) and second at qjjiy, —rj). 

Let T G C be the set of directions rj such that {y, —rj) is tangent to H at its points of 
intersection. 

U {Vi ~v) E C° = C\T, let t]j{y, rj) < t'jj{y, rj) be the times of intersection with H . 

The reason for the minus sign in —rj is that the transfer map follows the trajectory defined 
by {y^rj) in the backwards direction (see ( [T3l) . The time reversal is the map rj — )• —rj. Since 
H is smooth, the time functions t^{y,rj) are smooth functions of their arguments in C°. 

3.2. Formulae for the branches of the transfer map. In this section, we give explicit 
formulae for the branches of the transfer map away from T. They are needed in ^ in the 
proof of the genericity of the condition on H. 

We define branches Tq^ {y, rj) and Tq^ {y, rj) of Tq^ for {y, rj) ^ T corresponding to the 
two intersection points with H. For {y, rj) G C", we define the once if -broken trajectories 
EHj(t,y,rj) with breaks at q'^iy^v) ^^^ intersection times t-'j^{y,rj) by 

r E{t,y,-rj), t<t^{y,rj); 

EH,jit,y,rj) := < 

( E{t'H{y.v),y,V) + {t-t'^{y,rj))C'^{y,-rj), t>t^H{y^v) 
where 

C^{y, -v) = -ay, -v) + 2(C ■ ^,^(,,,))^,^, (,,,)• (38) 

As above, the —rj is needed since the trajectory proceeds backwards from the initial condi- 
tions. The unit normal is understood to be on the -|- side of N*H. 

After the break, the trajectory proceeds on the straight line E(t-'j^{y,rj),y,rj) + (t — 
'I^Hiy,v))CHiy,v) until it intersects the boundary again at a point Qnaniy^v)- We denote 
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by t-'fjQQ{y,v) the time at which this occurs. The following is a more detailed version of 
Definition |2j 

Definition 3.5. Define (see IfS^) ) (lHQfi{y,v) ^^ be the solution of the equation 

(iH{qH,dniy^ v)) = Eit^H^y^ v),y, v) + (^,9^(2/, v) - tuiy. vJKhiy^ v)- 

If(lH,dniy^v) ^ S, put 

^miy^ v) = WH,dniy^ v), T^T{q^^^^{y,v)fiqH{q^H,dniy^ ^))' ^i^niy^ v), y, -v)), 

Also define the correspondence 

^miy^v) = {^dn(y,v),^m{y,v)}- 

In the notation of [Z] , this double- valued correspondence would be written in the additive 
notation, 

^dniy^v) = r9^\y,v) + ^dn{y,v)- (39) 

If {y, rj) G T, then its trajectory runs tangent to H, and it proceeds on a straight line 
to dQ. We then define Tq^ (y,"!]) = '^ad (y^v) = PhJiV)) the usual billiard map. Thus, the 
graph r H-.j C B*dQ x B*dQ of the indicated partially defined symplectic map is given by 

r H =r H:l Ur H;2, (40) 

where the two graphs on the right intersect along their common values for (y, rj) G T. 

As in the discussion of /3, we note that Tq^ and Tq^ are piecewise smooth maps. It 
follows easily that the domain of definition of Tq^^ is C\EJj^-{Tj) and each map is smooth on 
the interior of its domain. 

In the case of more general hypersurfaces, billiard trajectories from dVt may intersect H in 
more points, and may intersect both tangentially and non-tangentially. As mentioned above, 
it is unnecessary to deal with such complications since the QER theorem may be proved for 
symbols a with arbitrarily small microsupport on B*H . 

3.3. Symbols of Nqq^ and Nqq . We first recall the definition of the symbol of a semi- 
classical Fourier integral operator F{X) e I^{M x M,C) where A > Aq is the inverse semi- 
classical parameter. It is a smooth section of the bundle fti/2 ® L of half- densities on C and 
the Maslov line bundle. 

In a local Fourier integral representation F(\){x, y) = (2nX)^"- J^„ e^*'*"^*^^'^'^)a(a;, y, ^, \)d^ 
with a(x, y, ^, A) ~a^oo J^TLo cij{x, y, ^)\~^ , aj{x, y, ^) G Sil, the symbol is transported from 
the critical set 

C^ = {{x,y,0 : d^v{x,y,0 = 0} C M x M x M" 
by the Lagrange immersion i^ 

i^:C^-^C, {x,y,^) ^ {x,ip'^,y,-ip'y). 

On Cip we form the Leray density 

J dxAdyAd^ 

^^ ~ d(|f ) 



22 JOHN A. TOTH AND STEVE ZELDITCH 

At a point (xq, ^o, Vo, Vo) ^ C, the symbol of F is defined as 



(^Fixo,^o,yo,r]o) = i^*ao^Jdc^ 



Lemma 3.6. The principal symbol of Nq^ is the half density on the graph of Tqq given at a 
point {{y,r]),Tg^{y,r])) G T^h with {y,r]) G B*dn by 



'an 



i{y,v)) V^'^i, , |i/2 



-(iVa"n(A))(y,r,)=(-P^) \dyd. 



li^aniy^v)' 

in the symplectic coordinates {y,ri) of B*dQ. 
Also, 

cr{\-'N^^^{X)){y,r^) = i{j{y,r^)MTUy,v)Y'\dydv\'^'. 

Proof. The calculation of the symbol of Nq^{X) is similar to that of Ng^{X) in Proposition 
6.1 of |HZ] . and we follow that calculation closely, referring there for some of the details 
and emphasizing only the steps where something new occurs. From the asymptotics of the 
Hankel function, 

oo 

i=o 

the principal symbol of Nq^{X) at the billiard Lagrangian is the same (up to an eighth root 
of unity) as that of 

where q{y) G dQ,qH{s) G H, and where diyy\q{y) — qH{s)\ is the directional derivative taken 
in M" X M". Hence the symbol is given by 

\qiy) - qnisr^'-'^^'d^Miy) - ^His)\\dyds\'^\ (41) 

where we use {y, s) G W^^^ x M"~-^ as coordinates on the graph of rj^. 

To express the graph half-density in terms of dydt] we need to express ds in terms of dr], 
keeping y fixed. Along the graph of Tq^, rji = dy.\q{y) — qni^)] where now the derivative is 
only along dfl. Hence, 

\dv\ = det (g^F-Q^ My) + y^ei) - {qn^s) + S'.e'JM \dsl 

where e, is an orthonormal basis for Tydfl, and e^ an orthonormal basis for TgH. This is 
the same type of determinant as in [HZ], Proposition 6.1, except that qnis) G H rather 
than in dfl. As in (loc. cit.) we first compute the determinant in the two dimensional 
case. We choose coordinates so that q{y) = (0,0), ^^(s) = (0,r), ei = (cos a, sin a) and 
62 = (cos/?, sin /3). Then the determinant is 

I (F cos a — S'cos/3,Fsin« — r — Ssmf3)\g^Y=o ' 
which equals 

r"^cosacos/3 = \q{y) - qHis)\'^d^ \q{y) - qHis)\d^Jq{y) - qnis)]. 
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Taking this number to the power —1/2 and substituting in ( HTl) . we find that the factors of 
Wiy) ~ Qnis)]'^^"^ cancel, and that there is a half power cancellation in d,yy\q{y) — qnis)], so 
that the symbol is given by 

^.Jg(y)-gH(^)h^/^|,^,^|V2 ^ ( ^%^Y' \dydv\'^' ■ (42) 



duMy)-iHis)\J Kiir^niy^v)) 

To complete the calculation, we observe that dujq{y) — qH{s)\ = 7(s,r), and diyy\q{y) — 
1h{s)\ = 'y{y,ri), and this gives the symbol of Nq^ in dimension two. 

In higher dimensions, as in |HZ] . we use the subspace Tq(^y)dfl fl Tq^(^s)H fl £^g where ig^s 
is the line joining q{y) and qnis)- We then introduce the same coordinates as in [HZ] and 
find that in the n-dimensional case, the determinant equals r"*^""*"^) cos a cos 7 where ei = 
(cos a, 0, sin a), 62 = (0,1,0) are the first two elements of an orthonormal basis of Tg(^y)dfl 
and e[ = (cos 7 cos /3, cos 7 sin /3, sin 7), e'a = (— sin/3, cos/3, 0) are the first two elements of 
an orthonormal basis of Tq^(^s)H. The factors of |g(y) — qH{.s)\ again cancel out due to the 
factor of r"""*"^. The cosine factors also produce half-cancellations as in the two-dimensional 
case, leaving the stated symbol. 

The calculation of the symbol of A'^^'" is similar except that the kernel has one further 
derivative, and the leading term arises by differentiating the exponent in v. So its symbol 
is that of A^^ multiplied by iXdy^\qH{s) — q{ii)\ where s E H,y G dfl and the derivative is 
computed in M" x R". In terms of the symplectic coordinates {y,f]) on the boundary, the 
symbol of Ng^ is the product of the symbol of A^|^ by the additional factor of 



i>^duMy) - Qnis)]. 



This completes the proof of Proposition 13.31 Q.E.D. 



D 



4. F(A) = N^^{XyOp,{a)N§,{X) and F'^(A) = iV|,-(A)*Op,(a)iVf^-(A) 
The purpose of this section is to prove that the compositions 

F(A) := N^^{XrOp,{a)N^^{X), resp. F-(A) := N^^^{XrOp,{a)N^^{X) (43) 



have the decomposition as described in fITTj) . The precise statement is given in Proposition 
14.21 The details and statements are the same for -F(A) and F'^{X) except for the calculation 
of the symbols. Hence, we sometimes give them only for F{X). 

By Proposition 13.31 and the calculus of Fourier integral operators, F{X) would be a Fourier 
integral operator associated to the composite canonical relation (r^)*or^ if the composition 
were non-degenerate (or clean); see [DSj . The composition is degenerate along the tangential 
set, and so in analogy with the boundary case in |HZj we must take care to avoid both the 
tangential and singular sets. We handle these problems by using appropriate cutoffs to 
remove these thin subsets. Thus, we first consider A-pseudodifferential operators Opx{a) G 
Opx{S°f{T*H X [0, Ao^])) with supp a f] T^^iTUi: U ^ ) = for A; = 1, 2; here as above, 
Tg^ : B*dQ — )■ B*H is the transfer billiard map. As in [HZj . a density argument shows 
that such operators suffice to prove Theorem H] in the general case. With this in mind, we 
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henceforth assume that the symbol a G 5*^/ {T*H x [0, Aq^]) satisfies the following support 
condition: For arbitrarily small fixed e > and any (sq, Tq) G rj^( TU ^ U B^dQ ), 

a(s,r;A)=0, (44) 

for all (s, r) G i?*-?/" with |s — sqP + |r — tqP < e"^ and A > Aq. 

4.1. The canonical relation of F{X). Temporarily ignoring the complications involving 
singular and tangential rays, the phase of the composition is \i/(s,|/,?/'); = \q{y) — (1h{s)\ — 
\Qh{s) — q{y')\ with qnis) G H,q{y),q{y') G dQ. Hence the critical set of the phase "^ is 

C* = {{s,y,y'); V.(|g(y) - qnis)] - \qH{s) - q{y')\) = 0} 

(45) 
= {{s.y.y'); {r{qH{s),qiy)),Ts) = {riqHis),q{y')),T,)}. 

The canonical relation of F{X) is then 

Tf{x) = {{y,Tryr{q{y),qH{s)),y',-7ry,r{q{y'),qH{s)); {s,y,y') G C^}. 

One component of Tp(^x) consists of the diagonal Aqxc C C xC (see Definition 13. 4p . Under 
the convexity assumption on H, this component occurs with multiplicity two since there are 
two points of H where the ray my intersect H and then refiect to its initial position. Two 
other 'components' of ^f{x) (which meet along rays tangential to H) come from once-broken 
trajectories from (y, //) G B*dQ to {y',^]') G dQ with the break at qnis) G H. There exist 
(at most) two possible break points corresponding to the (at most) two possible intersection 
points of a line with H. 

We denote the partial symplectic maps defined by the once-broken trajectories by (3}^, 
resp. Pjj, where the index indicates the first, resp. second point of intersection with H. As 
in the introduction (see flT^ we refer to these partial symplectic maps as transmission maps. 
Summing up, we have 

Lemma 4.1. Assuming H is convex, the canonical relation underlying NQ^{X)*Op{a)Ng^{X) 
(and ofN^^{XyOp{a)N^^{\)) equals 

which consists of three branches: 

• The identity branch on C (with multiplicity two); 

• The graphs of (5\j and Pjj; the index specifies the point of intersection with H . 

The case of general smooth if is of a similar nature but with more branches. We therefore 
omit the details. 

4.2. The decomposition of F{X). We now give a precise proof that deals with the A- 
pseudodifferential factor, the tangential degeneracy and the singular points. We continue to 
assume that H is convex, but again, this is for notational simplicity. The A-FIO contribution 
to N§^{XyOpx{a)N§^{X) is discussed in subsection gXH 

Let 6{e) > be small (to be specified later on) and let x ^ C'^(M) with xi^) = 1 when 
|x| < 1/2 and x{^) = for |x| > 1. We define the rescaled cutoff function Xs{e)(yX) := 
xisfs))- '^^^ point of the next Proposition is to show that by choosing symbols a satsifying 
(jH]) (ie. with support disjoint from the image under the transfer map of the generaltized 
grazing and singular sets), the corresponding -Fi(A; a, £:)-operator in the decomposition 
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is, modulo a residual operator, a A-pseudodifferential operator of order zero with principal 
symbol (Tj^)*ao ■ |o'(A^^)p. Following [HZJ we say that a A-dependent operator R is residual 
provided it is smoothing and \d'^d^R{x,y]\)\ = Oa^/3{X~°°)- The F2(A; a, £)-operator in the 
decomposition (l48l) is studied subsequently in subsection l4.2.1[ It is a sum of two zeroth-order 
A-Fourier integral operators. 

Proposition 4.2. Suppose that H is convex and that a e S^'°{T*H x (0, Aq"^]) with a ~ 
Yl'^o^j^'"' satisfies the support condition ( f^^P - Then, for any e > 

(zXjA)*OpA(aXn(A) = OpxiaoiT^niy,v)) " liv^v) ■ l'\r^niy,v))) 

+ F2i(A; a, e) + F22(A; a, 5) + R{X; e), 

{tz)X-'N^^^{XrOp,{a)N^^^{X)=Op,{ao{T^niy,v)) " liv^v) ■ lir^niy,v))) 

+ F^,iX;a,e) + F,',iX;a,e) + R''{X;e). (46) 

/n 6ot/i cases F2i(A;a, e) (^resj». F2''i(A; a,e)j and F22(A;a, e) (resp. F22{X;a,e)) are X- 
Fourier integral operators of order zero with K{F2i{X;a,e)) = K{F2i{X;a,e)) = /3^ and 
K{F22{X;a,e)) = K{F^2i>^;a,e)) = Pjj andmax{\\R{X;e)\\L2^L^,\\R''{X;e)\\L2^L^) =0,{X-^). 

Remark: The case of a general smooth H is similar but with a number of branches depending 
on {y,ri). We omit the details. 

Proof. Consider the operator A^^(A) : C°°((9f2) — )■ C'^{H) with Schwartz kernel N{qH, q; A) = 
du Golgn, q] A). By an integration by parts argument, it follows that 

wf;^{n^^{x)) c B*dn x b*h. (47) 

Using X5{e) we decompose the Schwartz kernel oi NQ^{X)*Opx{a)NgQ{X) into near-diagonal 
and complimentary terms: 

N^nW*Opxia)N^^{X)iq,q') = F^{X;a,e)iq,q') + F2{X;a,e)iq,q'); {q,q') E dn x dn, (48) 

where, 

Fi(A;a,£)(g,g') := iV|[,(A)*OpA(a)iVj'n(A)(g, g') x X5(e)i\q ~ q'\) (49) 

and 

F2(A;a,£)(g,gO := N^^{XrOpx{a)N^^{X){q,q') x (1 - ;^,(,))(|g - g'|). (50) 

We also introduce a new cutoff E^ G C^{Bl_^_^H) with S^d/, r^) = 1 for (y, 77) G Bl_^^,2{H). 
Then by the calculus of wave-front sets, we have 

F,{X;a,e){q,q') = N^^{Xr0px{E,)0px{a)N^^{X){q,q')xx5ie){\q-q\) + 0,{X~n 

uniformly for (g',g) G Sfi x 9i7. By fj47|) it suffices to assume that a G C^{Bl_^_^H) and 
satisfies the support condition (jUj) and we will do so from now on. 

From Proposition 13.31 we know that Nqq{X) is a zeroth-order A-Fourier integral operator 
with canonical relation T^h . Differentiation of the phase function "^(yjy^s) = Iqnis) — 
q{y)\ — |i?h(s) —q{y')\ in s and a repeated integration by parts implies that for some constant 
C>0, 

iyF{(Fi(A;a,£)) C {{y,^,y',v) e B*dQ x B*dn; \y - y'\ < 6{e), \^ - v\ < C6{e)}. 
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Since (see section KA}i . 

{Tp.^uTp2^)nAcy,ccrug, 

and dist(supp a, Tq^{T, UTuQ)) > e, it follows that by choosing 6{e) > sufficiently small, 



^H I^H' 



Thus, in view of Lemma 14.11 

WF'^{Fi{\;a,e))cAco,,co. (51) 

By the semiclassical Egorov theorem |DS] . Fi{\;e) G \E'°^(9n), and moreover, 

Fi(A;a,e) = Op,( (rfj*ao ■ |a«^(A))nn + ^(A-^)^.^^.. (52) 

From Lemma ESI |c(iV'g^(A)(y,?7)p = ■j{y,ri) x 'J'^ {TQ^{y , r])) and so the formula in (1521) 
implies that 

F,{X;a,e)=Opx{ao{T^^{y,v))-l{y,v)-l~\r^niy,v))) + OeiX-')L^^L^. (53) 

In case (ii) one makes the analogous decomposition 
N^n'{XrOp,{a)N^^^{X){q,q')FaX;a,e){q,q') + F,^{X;a,e){q,q'y, {q,q') e dQ x dn (54) 

where, 

F^{X-a,e){q,q') := N^^^{XrOpx{a)N^^^{X){q,q') X Xsie){\q - q'\) (55) 

and 

F^{X;a,e)iq,q') := iVj[r(A)*Op,(a)iVj[f (A)(g, g') x [I - ^sie))i\q - q'D- (56) 

The same reasoning as in case (i) implies that 

F^{X;a,e)=Op,{{T^^rao-\a{N^^^{X))n + 0,{X-')L^^L^ (57) 

and by the symbol computation in Lemma 13.61 one gets 

Fr(A; a, e) = Op^{ a^{T^^{y, r/)) ■ 7(ri^^(y, r/)) ■ 7(1/, r^) ) + 0,{X-^)l2^l^. (58) 

4.2.1. Analysis of the F2{X; e)- operator. In this section, we prove: 

Lemma 4.3. F2i{X;a,e) (resp. F22{X;a,e)) are X-FIO's of order zero with K{F2i{X;a,e))l3}j 
(resp. K{F22{X]a,e)) = (3^). Their principal symbol (on the respective branch of the corre- 
spondence) is given by 

dvy^\q{yi) -qH{s)\\^''^ (d^^^\q{y2) -qni-sy^ ^'"^ 



"*' d.Mivi) - qH{s)\ ) \d,M{y2)-qH{s)\) \^=^H-Ay^m) 

-air-,iy2,V2))( /^f/^^^ j^^V /^f/^^^ )^^^; iy.,v.) ^ ^UvuV.), k ^ 1,2. (59) 
^7(^ah(z/i,^i)^ ^l{ron{y2,V2)^ 

Here, s = qH,k{y, y') & H;k = 1,2 is the point on H joining y, y' G dVt under the H -broken 
billiard map /3|^ :; A; = 1, 2, (see section\^. 
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Proof. Because of the cutoff (1 — X5(£))(|5'(l/)— Q'd/OI) i^ ^^e amplitude of the F2{\; a,5)-kernel 
the canonical relation corresponding to -F2(A; a, e) does not intersect the diagonal AB*dnxB*dn 
and so, by the argument in section [3] it must consist of the reflection canonical relations Pjj 
(resp. f3jj ). The corresponding quantizations are the A- Fourier integral operators -F2i(A; a, e) 
(resp. -F22(A; a, e)). Modulo pointwise uniform 0(;(A"°°)-errors, the Schwartz kernel is locally 
given by the formula 

F2k{X;a,e){y,y') = {2nXr-' [ e^^[l''«(^)-''(^)l-l''-('^)-^(^')l] c(i/,i/', .; A) (1 - X5(e))(|2/ - 2/1) 

Jh 

X Xeiqnis) - QhAv^ y'))d(TH{,s) (60) 
where, c ~ Y^T=qCj\~^ with 

Co{y,y',s) =ao{qH,j{s),dsqH,j{s)r{qH,j{s),q{y'))) x {uy,r{qH{s),q{y))) 

X {J^y', r{qH{s), q{y'))) ■ bo{q{y), qH,j{s)) ■ bo{qH,j{s), q{y')) (61) 

n— 1 

where, &o(<?h(s), g(2/)) = \qH{s) - q{y)\~^ ■ 

The symbol is computed by taking the product of the principal symbols of the factors 
given in Lemma [3.61 

D 

Combining f l52|) with Lemma [4.31 completes the proof of Proposition 14.21 in the case of the 
N^nW*Op,{a)N^^{X)- operator. 

The following subsection completes the proof of the Proposition. 

4.2.2. Modification for N^^{\)*Op{a)Nf^^{X). 

Lemma 4.4. F2i{X;a,e) (resp. F22{X;a,e)) are X-Fourier integral operators of order one 
with K{F2i{X;a,e)) = (3^^ (resp. k{F22{X] a, e)) = (3q^ ). Its symbol 

- {duyJqiVi) - qHis)\d^Myi) - QHis)\Y^'^x{d^Jq{y2) - qH{s)\d^Jq{y2) - qH{s)\Y^^ \s=gi,.^^(y,,y,) 
is that of F2J multiplied by 

-(iXd^skiyi) - qH{s)\){-iXduMy2) - qH{s)\)\s=g„._^(y,,y^). 

In symplectic coordinates, it is given by 

- 7(2/i,^i))7(^i^(2/i,^i) 7(2/2,^2))7(^^(2/2,^2) ; (2/2,?/2) = /3h(2/i.^i). ^ = 1.2. 



The proof is the same as for Lemma 14.31 and is therefore omitted. This completes the 
proof of the Proposition. D 

5. Proof of Theorem [U 

It is important to observe the sign difference between the symbol of the FIO part of 
A'J^(A)*Op(a)A^K^(A) and its i/-partner. In the calculation of the first symbol, the normal 
derivative with respect to s & H comes about when we change from (s, s') coordinates to 
(s, r) coordinates. So this normal is symmetric with respect to the two sides of H (at each 
intersection point). But in the Nq^" operator it is assymetrical. This explains why the F2 
operators have opposite signs in the two A^*A^ terms. 
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We note also that in the case of the A-pseudo differential F^{X; a, e)-term, yi = y2 and so, 
yi and 1/2 trivially lie on the same side of H relative to the normal vector J^s=qH.,,{yi,y2)- On 
the other hand, in the A-Fourier integral case of -F2^(A; a,e), yi and 2/2 he on opposite sides 
of H relative to the normal vector Vs=qHk{yi,y2)- This creates the extra minus sign in Lemma 
14.41 which allows for the cancellation between the (appropriately weighted) F2 and Fg^'-terms. 

We use these observations to reduce the proof of Theorem [1] to that of |HZ] : 

Lemma 5.1. The principal symbol of the Fourier integral operator part of 

<n(A)*OpA((l - |rnaX^(A) + \^'N^^{\yOp{a)Ni^{\) 

equals zero. 

Proof. The principal symbols of the F2{X) terms of the two operators are calculated respec- 
tively in Lemmas 14.31 and 14.41 As discussed in the proofs, the second (Neumann) sym- 
bol equals the first (Dirichlet) multiplied by the factor of du^\q{yi) — qH{s)\)dyJ\q{y2) — 
to('S)|)|s=g^.fe(j/i,y2)- By definition, 5(1/2) — qui^s) and qiyi) — qnis) have the same tangential 
projection to H. But they lie on opposite sides of H and have opposite normal projections. 
So this factor equals — cos'd{yi, s)^ where ■& is the angle to the normal. In symplectic coordi- 
nates of the projection (s, r) on B*H, cos'i9(?/i, s)^ = 1 — |rp. ) Hence, if we take a diagonal 
pseudo-differential operator A2 with o"Aii(s,r) = (1 — |rp)a(s,r) and (7^22 (■S) t) = a{s,T) for 
all [s, r) G B*H, and if we normalize the second term by dividing the Neumann data by A, 
then the Fourier integral parts of the NQ^{X)*Opx{a)Ng^{X) compositions cancel each other. 
Then we are left with only the pseudo-differential parts of each composition, completing the 
proof of the Lemma. D 

Thus, by Lemma 15.11 and Proposition 14. 2[ 

(An(A)Mf ,Mf )i2(^) + X-^{A22i\)u^,u^)LHH) 

= 2 {Op,{ {T§,ra ■ {ri^r^ ■ ^)ulu\)L.^en) + 0{X-')\\u\\\l,^,^^ (62) 

and the statement of Theorem [1] then follows immediately from ( 162|) and an application of 
the boundary quantum ergodicity result of |HZ] . D 

6. Boundary Weyl law: Proof of Theorem [9] 

We now prove a key result underlying Theorem H] on quantum ergodic restriction for 
Dirichlet data. The proof hinges on an averaging argument in ^ and on a boundary local 
Weyl law. The purpose of this and the next section is to prove a weak version of a boundary 
Weyl law for A-Fourier integral operators. The presence of the boundary causes a number 
of technical complications, but the result seems to us of independent interest. This section 
is independent of the choice of H. 

To prepare for the boundary result, let us first recall the case of manifolds without bound- 
ary and eigenfunctions of a Laplacian [Z] . The Weyl law for classical Fourier integral oper- 
ators F : C°°{M) -^ C~(M) of order zero says that 

J2{F^x,,^x,)-^ f (T^{F)d^i, (63) 

.\ ^ \ 'J rpflA'T-'* AT 



^(A) ^,— ^ Jr^nA^: 
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where the integral is over the intersection of the canonical relation Tp of F with the diagonal 
of T*M X T*M, i.e. the fixed-point set of the symplectic correspondence underlying F. The 
right side is zero unless the fixed-point set has dimension m = dim M, the leading order trace 
sifts out the 'pseudo-differential part' of F. The fixed point set is assumed to be almost-clean 
in the sense defined above Theorem [H 

The purpose of this section is to prove an analogous result for semi-classical Fourier integral 
operators acting on boundary values of eigenf unctions. Our goal is to prove (Theorem [9]) 
that 

^^^(A)(A) := Yl (i^(A,)<.,<.) = 0(1) (64) 

j:Xj<X 

if |Sj7(A)| = or if the quantitative almost nowhere commuting condition holds. 

The study of these A-FIO Weyl sums introduces two new aspects to the local Weyl law. 
First is the semi-classical aspect, which would also arise in local Weyl laws for matrix elements 
{F{Xj){pxj,ip\^) of semi-classical FIO's relative to Laplace eigenfunctions. As is usual with 
semi-classical FIO's, we 'homogenize' the traces by taking the Fourier transform in the 
semiclassical parameter A and integrating over the dual frequency variable. 

The second novelty is the restriction to the boundary. The semi-classical parameter is 
the interior eigenvalue, not the eigenvalue of an operator on the boundary. The restriction 
operator to the boundary resembles a Fourier integral operator, except for the complications 
of grazing rays, corners and so on. These play a minor role in the trace since they occur on 
sets of measure zero and are handled as in |HZj by introducing appropriate cutoff operators. 
We omit some of these technical details for the sake of brevity, since they are the same as 
in |HZ] . In the subsequent article |TZ3] . we prove a more complete boundary Weyl law for 
A-FIO 's without assuming measure zero conditions. 

Before stating our result, we introduce some further notation and assumptions. We assume 
that there is an coveing IJ^i ? (^) — ^^ with U C M"^^ open, such that in terms of the 
local parametrizations y t— )■ q^~'\y);j = 1, ..., K, the phase function has the form 

Hq^'\y),q^^\y'),e) = {y,e)-^^^\y\e), ^^^^ G C^{U x W'^). (65) 

The associated canonical relation (135|) is the graph of a symplectic correspondence 

kf : B*dVL -^ B*dn. (66) 

This is sufficient for our purposes since the semi-classical FIO's of interest here are the 
pieces of F{\)N§^{XyOpx{a)N^^{\). We denote by r^o : S*n9n -^ S*nan the reflection of 
inward pointing covectors at the boundary with respect to the interior unit normal given by 
i"oq(C) = C + 2(C,^'g)^'q; q £ dQ. Also, the map tit '■ Sq^Q — )■ B*dQ denotes the canonical 
tangential projection along dfl and we recall that C^ denotes the broken geodesic (billiard) 
flow (129|) . From now on, we drop the (j)-superscripts in the phase functions ( 165|) with the 
understanding that computations are local. 

6.1. Proof of Theorem [91. 

Proof. In the following, we assume that dQ is strictly convex. We treat the general non- 
convex case in subsection 16.71 
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Following the strategy of the Fourier Tauberian theorem |SV] , we calculate the principal 
term of the Weyl asymptotics of order A" by computing the singularity at t = of order 
(t + iO)-" of 

l\t):=J2e'''HFiX,)u\^,u%). (67) 

j 
We note that by adding a constant multiple of the identity operator, we may assume that 
the Weyl sums are monotonically increasing, hence the Fourier Tauberian theorem applies 
(see [ZJ for the parallel statement in the boundaryless case). 
To calculate the principal singularity of ( |67|) at tO, we write 



F(A) = ^ f e^'^F{s)ds, 

^T^ Jr 



where 



Jr 
Then F{s) is a homogeneous Fourier integral operator with Schwartz kernel, 

F{s){q,q') = {2nfir-' [ [ e^'^(*(«'«''^)-)A(g, g', ^,/.) d^rf/i, (68) 

Jr+ JR"-i 

where A{q,q',e, fx) ~ EJlo^i(?' ^'' ^)/""^' ^i^h Aj E C^{dn x dn x M"-i). We then have, 
5^e^*^^(F(A,)n^^.,n^; = / 5^e'^^(^+*)(F(.)u^^,<,)d., 

or equivalently, 

I\t) := Tr I ^l^F{shenUit + s) ds. (69) 

Jr 
The presence of the boundary restriction operator 7^^ has the effect of restricting kernels 
to the boundary. To see this, we note that if jh denotes the restriction operator to H, 
IhI = f\H^ then Ynf = f^H^ since {Yuf^d) Jh fd^^- Hence, in the Neumann case, writing 
U^{s) = 'ydnU{s)'ysn', s G M for the boundary restriction of the operator U{s), 

I\t)= f f [ U\t + s,q,q')F{s,q',q)dsda{q')da{q). (70) 

Jr Jan Jan 

In the Dirichlet case, restriction first takes normal derivatives. In the following, we wish to 

compute the coefficient of the large (t + iO) ""-singularity of I^{t) as t — )■ 0. 

6.2. Microlocal cutoffs. To simplify the trace we introduce some microlocal cutoffs: Since 
both the singular set, S, and the grazing set, Q, are closed and of measure zero, by the 
C°° Urysohn lemma, for e > arbitrarily small, we let Xs ^ C^{B*dVt) be supported an 
^-neighbourhood of E U ^.. Such cutoffs are necessary but quite standard [HZ]. 

In addition, we will also insert a cutoff in the integration in s in 07Uj) near s ~ $ 
and localize the trace in time. This cutoff is more novel and is also important since it 
allows us to ultimately replace U^{t + s) by the corresponding finite-time reflection (ie. 
Chazarain) parametrix in the formula fITOl) for the boundary trace, I^{t). We define the 
A-microlocally cutoff operators -Fe(A) := Op\{l — Xs)F{X)Op\{l — Xe) and U^it + s) : = 
Op\{l — Xe)U^it + s)Opx{l — Xe)- By standard composition calculus, one can write the 
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Schwartz kernel F,{\){q{y),q{y')) = {27r\)''-^ J^^^, 6'^'^'^'^^^^'''''^'^'^^ A,{q{y),q{y'), 9) dO where 
Aeiq{y), q{y'). e) = Aiqiy), q{y'), 9) ■ {I - XeMv), dy^) ■ (1 - XsMv'), dy,^) + 0{\-'). The 
precise statement of the A-microlocahzation of the boundary trace is given in the following 

Lemma 6.1. Let x e C^(ffi) with x{s) = 1 near s = and for any e > let Xe ^ C^{B*dn) 
be supported in an e -neighbourhood of B^dQ U Q. Then, 

I\t)= Iff f f /7,^(t + .,g,g')e^'^W''^'^)-^)A,(g',g,e,^) 
JR JR jR"-i J an Jan 

X xi-s - $(g, q', 9)) {2TTij,)''-^da{q')da{q)d9dsdfi + R{t) + R^it), 

where R{t) G C°°(Ri) and R,{t) = 0{e)t-'' for \t\ small. 

Proof. From L^-boundedness of F{X) and the boundary Weyl law (see [HZj Lemmas 7.1, 9.2, 
Appendix 12 as well as Lemma [8.51 in section [8^ of this paper), it follows that 



iV(A) 



Y,{OpxMe)n>^^)OpxMeK^<,) 



\i<\ 



< -£-J2{0px{Xsr0px{Xe)n%,u%) = 0{e). 



NiX) 



A,<A 



Similarily, by applying Cauchy-Schwartz and an argument like the one above it follows that 
the composite operators Opx^{xe)F{Xj) Opx^ {l - Xe) and Opx^{l - Xe)F{Xj)OpXj{Xe) each 
have 0{e) boundary traces (see also [Z| IMP] ). 

Thus, substitution of the the integral formula (l68l) for -F(s, g, q') in the formula for I^{t) 
gives 

I\i) = L Ian Ian U\t + s, g, q')F{s, g', q)dsda{q')da{q) 

= /m+ /r /m"-i Ian Ian U'eH + s, ?, g')e^^^*^'^'''^''^"^^^(?', ?, ^, ^^) (2^"-^ da{q')daiq)d9dsdfi + R,{t). 

We now cutoff in the s-time variable. First, note that since A^ is symbolic in /i, 

9;A,(g,g',^,/^) = 0,,,(/x-") 

uniformly in the (g, g', 5)-variables. 
We make the following decomposition 

I\t)= f f f f / f/,^(t + .,g,g')e^'^(*(«'''''^)-^M,(g',g,^,/i) 

Jr+ Jr jr"-! Jan Jan 

X x(s - $(g, g', 9)) {2nfi)''-^da{q')da{q)d9dsdfx + R{t) + R,{t), (71) 
where, j^ e^t'^Re{t)dt = C(e)/i"-i and 

^W=/ // / f U'At + s,q,q')e'>^^''^^''''''^-^^A{q',q,9,fi) 
Jr+ Jr Jr"-i Jan Jan 

X (1 - x)(s - <l>(g, g , 9)) {2'K^iY-^da{q')da{q)d9dsd^i. (72) 
By integration by parts in /i and the fact that d'^A = Oaifi'"), it follows that for any A^ > 0, 

m= [ II I / f/.'(t + 3,g,g')e^^(*(^''^'')-^M,(g',g,e,/i) 



e 

"-1 Jan Jan 



32 JOHN A. TOTH AND STEVE ZELDITCH 

X (1 - x)(s - $(g, g', e)) X Cjv(|sr^/i"-^"^) da{q')da{q)dedsd^i. 

Since the integrand is absolutely convergent, one can take the t — ?■ limit inside the integral 
and this shows that R{t) is continuous near t = 0. Continuity for the derivatives of R follows 
a similar fashion: Since dtU^{t + s) = dsU^{t + s), one integrates by parts in s. Each 
s-dertivative of the exponential e*^*^*"*-* creates a power of fx. However, this is killed by 
the (9Ar(|s|^/i'^~^)-factor since A^ > is arbitrary. The result is an absolutely convergent 
integral. By iterating this argument, one can differentiate to any order in t inside the 
integral and the result is absolutely convergent. It follows that R{t) G C°° for |t| small and 
this finishes the proof of the lemma. D 

Thus, it is enough to consider from now on the principal part of the trace, I^(t), which by 
Lemma [6. H is given by 

X x(s - $(g, g', 9)) {2'K^)'^~^da{q')da{q)dedsd^. (73) 



sing \ J ■ I I II I e 

~ ' ~ I do. Jan 



The remainder of the proof consists of carrying out a singularity analysis of the integral 
for Igingit) in (173|) near t = 0. It is useful to note that since it suffices here to choose |t| < Eq 
with Eq > fixed arbitrarily small, \t + s\ < ^o + sup |$(g, q', 6*)! < oo for (g, q', 6) G supp A^ 
since supp A^ is compact. Thus, an immediate consequence of Lemma EH] is the identity 

/^(t)=/L,(t)+i?e(t)+i?(t), 

(74) 
JL,(t) = Tr j] U^{t + s) {T,^sFs){s) ds. 

Here, we recall that 

F,(g, q'- /x) := (27r/x)"~i f e*'^*(^'''''^)A,(g, g', O; /x) x{s - $(g, g', 9)) dO (75) 

and for e > arbitrarily small, Re{t) = 0{E)t-'', R{t) E C°°(Mt) for |t| small and / : = 

[-||$||loo -Eq, ||$||l«= +£:o]- 

6.3. Chazarain parametrix. Since U^{t+s) is A-microlocalized away from the set {B*dflx 

[Q U B^dVt]) U ([^ U B^dVt] x B*dfl), this allows us to substitute the boundary trace of the 
Chazarain parametrix for f/?(t + s) in the formula for /^j„c,(t) in (175]) . Here, we use the 
crucial fact that by Lemma 16.11 (see also fl7^ ). one need only consider finite times with 
\s + t\ < sup |$| +£0- This parametrix construction is well-known [GMtlPS] . However, as far 
as we know, the properties vis-a-vis boundary restrction are not available in the literature. 
Thus, we review here the aspects that are most relevant to the analysis of the boundary 
traces. 

It suffices to assume that dQ is C~ since WF;^{U^{t + s) o F,{s)) n {[B^dQ ug]x B*dQ U 
B*dQ X [B^dQ U Q]) = 0. Let dQ = {x eW' : f{x) = with df{x) ^ ior x e dft. Given 
{q,uj) e Slidn), we let < ti{q,uj) < t2(g,w) < t3iq,uj) < ... with tj G C°°{Sl{dn)) be the 
ordered intersection points of the bicharacteristic starting at (g,w). Thus, the tj;j = 1,2.... 
satisfy the defining equation 

/(t,(g,a;)) = 0;j = l,2,3,... 
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The remainder of the proof of Theorem [9] consists of substituting the boundary trace of 
the exphcit Chazarain parametrix (which we describe below) in (17^ and then carrying out 
a singularity analysis of I^{t) near t = 0. 

To describe the parametrix boundary trace, we decompose S^dQ = U^r^ where the F^ are 

sufficiently small open sets and let tj = min(g,j)grfe tj{q,u) and TJ = max(q_^)grs. tj{q,u). 
We choose T^ small enough so that 

^(fc) < j^ik) ^ ^{k) < yW < . . . . ^ _ ^^ 2, 3, ..., N. (76) 

Let Xk G C°°{S*dQ);k = 1,...,N, be a partition of unity subordinate to the covering by 
rk;k = l,..,N. 

One constructs a finite-time microlocal parametrix for the wave operator, U{s), on the sets 
Ffc; k = 1, .., A^. To describe the construction, consider first the case where s G (—260, A —Eq) 
where t[ := min(g,^)grj. ti{q,uj) and ^o > is a fixed small constant. We let 

[/f '+(s, X, y; A) := f/o o Op,(xfc)(s, X, y; A) = (27rA)" / e'^'^^^'^'^'^^sf (s, y, ^, X)d^, 

where Uq is the free-space solution of the wave equation with ip{s, x, y, ^) := {x — y,0~ ^\C\- 
Let X* {x) G M" — fi be the geodesic refiection of x G f2 in the boundary dQ and put 

f/f '-(s,x,i/;A) := UooOp,ixk)is,x*ix),y;X) = (27rA)" / e^^'^^^'^'W'^'^^sf (s,2/,e, A)^^ 



and for (s, x, y) G (— ^o, ti + ^o) x ^ x ^, define 

2t/f ) (s, X, y; A) := ?7o^'^'+(s, x, y; A) + f/o^'^'" (s, x.y; A). (77) 

Clearly, for all {s,x,y) G (— 2£:o,ii — £0) x ^ x ^, U{ ' solves the equation ^dsU[ ' + 

V^A^f/} = and since x* G M" — fi, by integration by parts in ^ it follows that for 
{x,y) eQ X Q, 

ui''\0,x,y;X) = OpxiXk){x,y) + R{x,y;X), (78) 

where R{x,y;X) G C°^ with \d^dlR{x,y; X)\ = 0„,^.(A-°^). Following pZ], we call such 

kernels residual. Simiarily, since d,y^[(p{s,x*,y,^) + f{s,x,y,^)]\x=x*=qedn = 0, f/| satisfies 
the approximate Neumann boundary condition 

d,^ui''\s,q,y; A) = R^is,q,y; X), {s,q,y) G (-£o,tf^ +^o) x 9f] x fi (79) 

with Ri{s,q,y; A) = (9(A~°°) residual. Also, it is clear that 

[Ui'^]\s, q, q'- A) = [U^'^^^]\s, q, q'- A); (g, g') e dQ x dQ. 

We microlocally cutoff f/| away from grazing and singular directions and define the mi- 
crolocal parametrix 

Ui'^is; A) := (Id - OpxiXe)) o u[''\s; A) o {Id - Op.ixe))- (80) 

The analogous, microlocally-cutoff wave operator is 

[/f )(s; A) := {Id - OpxiXs)) o U{s) o OpA(Xfe) o (Jrf - OpxiXe))- (81) 
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From (!78l) and (!79l) . both Ue (s; A) and f/f^ (s; A) are microlocal parametrices for the wave 
operator on supp Xfc(l — Xe)- Thus, t/j ^ (s; A) — Ue (s; A) is residual for s E (— 2£:o, 4 ~ ^o) 
and in particular, for (t + s, q, q') G (— 2£:o, n £^o) x dil x Sfi, 

[f/f ^l^lt + s, g, g'; A) = [f/J5]^(t + s, g, g'; A) + i?(t + s, g, g'; A) (82) 

where, R{-, -, ■; A) G C~( (-2eo, tf^ - Eq) x dn x dn) ) is residual for A > Aq. 

To construct parametrices up to finite, multiple reflections at the boundary, we choose 
M = M{k, e) large enough so that the finite time-interval, I of integration in (!74|) satisfies 

I C U,1i(f-\ - 2^0, f) - .o) [J-U-if^ - 2eo, -tf_\ - .o) with i'^ = and we let xg e 
C^(]R); j = 1, ..., M be a partition of unity subordinate to this covering. In view of fl52]) the 

parametrix for s + t G (— 2£:o7 A — Eq) is U{ J{s + t). For times t + sE{t[— 25o7 4 ~ ^o); 
one puts 

U^';^is + t,x,q;\) = J^ui';^is + t-T,x,z;X)-Ui';^iT,z,y;\)dz, r := \z - q{y)\. (83) 

Clearly, ^dsU2^ + A/A7f/2£ = ^^^ by (1821) . f/2£ (s + t) satisfies the Neumann boundary 
condition. 

Moreover, iteration of this construction gives the formula 

Uj'^e (-5 + t, X, y; A) = /^ f/}J {s + t-T,x, z; A) ■ f/-!\,e(r, z, y; A) dz, 

(84) 

s + tG(tf_\-2£o,tf -£o). 
Just as in the j = 1 case, it follows that for any j = 1, ..., M, and s + t G (4_\ — 2^0, ^ "^o)? 

[t/jj]^(s + t, g, g'; A) = f/f (s + t, g, g'; A) + i?,(g, g'; A); (g, g') G 9^^, (85) 

with Rj residual. We summarize the parametrix approximation of the boundary trace in the 
following 

Lemma 6.2. Given xf^ e C^{R);j = ±1,...,±M as above, it follows from (7^ and [M) 
that for \t\ sufficiently small, Igingit) = J2j ki-^j t(^) where, 

[/fY(t) .,= Trj]xf\s + t) {[U^'^]\s + t)oFe{s)) ds 

= Tr j\ xf\s + t) ([f/jj]''(. + t) o F,(s)) ds, 

where Uj J is the microlocal parametrix with Schwartz kernel defined in [84\ ) 

To carry out our analysis of the principal singularity of I^ingif) we use the following 
generalized stationary phase lemma ([SVJ Proposition 4.1.16): 

Lemma 6.3. Given ^(x) G C°°(R") with critical set Grit{^) := {x E W : V^^(x) = 0} 
anda{x) eC^(R'^), 

e^^*(^)a(x)rfx = f e^^*(^)a(x)rfx + o(l). (86) 
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The computation of the principal singularity of l''{t) amounts to a detailed stationary 
phase analysis of the explicit boundary trace of the A-microlocalized Chazarain parametrix 
given in Lemma [6.21 combined with multiple applications of ( 186|) . The analysis is somewhat 
lengthy since there are many critical points of the phases in the parametrix coming from 
various multilink bicharacteristics. We thus carry out the analysis systematically for the 
various links indexed by the the number of boundary reflections. 

6.4. Single-link contribution to /^^ (t). We start here with the computation of [/i]^(t) 

for small t with t G (— 260, niiufc t[ ]. This is the part of the trace that comes from single-link 
contributions. These are rays that are either trivial (ie. ones that intersect the boundary at 
only the initial point) or ones that intersect twice, with the initial and terminal points lying 
on the boundary. The main result of this section is 

Lemma 6.4. Under the assumptions in Theorem\^ the single-link contribution to the bound- 
ary wave trace is 

Ilit) = o{t--) 

for \t\ sufficiently small. 

In view of Lemma [6.21 substitution of the restricted parametrix [t/igX in [/} ]^{t) com- 
bined with the polar variables decomposition ^ = ru where (r, cj) G [0, oo) x S"^^ and the 
semiclassical rescaling ^ >-)• /x^, gives the formula 

[lf)]\t)= f [[ [ [ [ e'^[^^(^+*''?'^''")+*(«'''?'^)-^U,(g',g,^,/i)5f)(g,g',rw,/i) 
Jr+ J I Jm"-i Jan Jdn Js^-^ 

X x(s - $(g', q, 9))xir - 1) r""^ drdco Xi(s + t)deda{q')da{q)ds (27r/i)2"-i dfx + R{t). (87) 

Here, R{t) G C°^(M) for |t| small and the cutoff x(r — 1) is inserted in fl57|) modulo smooth 
error by integration by parts in s using the eikonal equation ds(f = r-\-0{t) with |t| ^ 1. Also, 
the integration in u is over the hemisphere 5*""^ since iyF{(f4 (s + 1)) C Sqq _^il x Sqq _^Q, 

for all s + t G /. From fl87j) it is clear that [JJ ]^(t) is a semi-classical (in /i) oscillatory 
integral with phase function 

^(g, q', r, w, 6, s; t) = np{t + s, g, g', u) + $(g', g, 6) - s. 

We apply stationary phase to (187|) . Since the integrand in (IHTj) is now compactly supported 
in (s,r), one can apply stationary phase in those variables first. Since 

ip(t + s, g, g, uj) = {q- g', w) + (s + t), 

the critical point equations are 

dr^ = ip{t + s, g, g', u) = 0, 

ds^ = rds(p{t + s, g, g', w) — 1 = 0. 
The non-degeneracy of the phase is clear from the identity 

drds'^iq, q', r, u, 6, s; t) = dsip{t + s, g, g', u) = 1 
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which follows from the eikonal equation. The critical points are 

r = 1, 

s{q,q',u) = -{q-q',uj) - t. 
The result is that 
[li^'m = L, /m^- Ian Ian V^ e»'^*e^'^[*(^'-.^)+<-^'-)]i?f (-(g - q', uj), q, q', uj) 

xMq, q', e, /i) x(-(g - g', oo) - $(?', g, e)) Xi(-(g - g', oo)) dujdeda{q')da{q) fi^^-^ d/i + R{t) 

We put ci := ^^ c[ and then from the last line of fl88|) it follows that 

JiR"-i Jan J an Js"~^ 

(89) 
where, A^{q,q',9; fi) ~ X]m=o ^e:"i('?' '?' ^)/^ "^ ^^"^ -^('^) ^^ °^^y singular of order 0(t "■^^). 
Here, we have also used that Bq = 1. 

We further decompose ci (/i) into diagonal and off-diagonal parts by introducing a further 
microlocal cutoff. We write 



ci(/i)=/ / / / ^e^'^[*(«'''''^)+<''-^'->U,,o(g,g',^)x(^-'|g-g'l) 

X Xi{~{(l ~ q' j^)) dcod6da{q')d(T{q) 



jR"-i Jan Jan Js"-'^ 

X Xi(^(Q' ~ q' , oo)) dujd9da{q')da{q, dcjd9da{q')da{q) 

=:cf(/i) + cf(/i). (90) 

6.4.1. Diagonal term. The term c^{iJi) sifts out the A-pseudodifferential piece of -F(A). In- 
deed, since |g — g'| < £, by Taylor expansion, 

{q{y)-q{y')^u) = {Ty>,u){y-y') + 0{\y-yf), 

and so, by making first the change of variables r : S^_^^ — ?■ 5""^, given by r(a;) = w — 
{vyi,uj)vyi =: r/ followed by another of the form rj i-t- rj{l + 0(\y — y'\)), it follows that 

cf(/i)=/ / / / e^^[^'^-'^(^'^)+<^-^''''>U,,o(g(y),g(2/'),^)(i + (^(l2/-y'l)) (91) 
jR"-i Jan Jan Jb"-i 

X X{£'^\y - y'\) l{ri)-^dr]dedo{y')do{y). 



QUANTUM ERGODIC RESTRICTION 37 



Here, we recall that 7(77) = yl — |?7P. An application of stationary phase in {y',0) implies 
that the RHS in flOTl) simplifies to 

(27r/i)-("-i) / / e^'^[<^'^>-^(^''')U,,o(g(2/),g(V,^),r7) {l + Oi\y-V,4j{y,v)\))l{vr' dyd^. 

(92) 
The assumption |Sp'(A)| = in Theorem M implies in particular that jF^^ fl A| =0. 
Consequently, the generalized stationary phase estimate ( l86l) applied to the integral in ( l92l) 
implies that for the diagonal piece, cf (/i), of the single-link contribution to /^(t), 

cf (/i) = o(/i-("-i)). (93) 

Thus, it follows that the contribution of the diagonal to the trace /^ingl^) i^ Io° e~**^cf'(/i)/i^*^"'~^^(i/i 
o{t~^) for |t|-small. 

6.4.2. Computation of cf'^^fi): This piece sifts-out the contribution to Ci(yu) of single links 
with initial and terminal points on the boundary. To compute the asymptotics of cf'^(/i), we 
apply stationary phase in to EE '8>'^~^. The solution of the critical point equation 

d^'^ = d^{q-q',u) = 

is the non-degenerate critical point uj{q,q') = —f-^\ (we note here that the integrand in 
Ci(/i) is supported on the set where \q — q'\ > s > 0). There is only one critical point, since 
the integration is only over u E S"~^ corresponding to directions pointing into the domain 
Q. The result is that 

cf (/.) = (27r/.)-("-i)/2 / / / e^^[*(''(^')''^(^)'^)-l''(^)-''(^')IU,,o(g(y),g(y)',^)(l + 0(/.-^)) 

Jm"-i Jan Jan 

x(l - xKe-'\y - y'\) xi{\y ~ y'\)da{y')da{y)d9, 
where, <^{q{y'),q{y),9) = {y',9) - ij{y,9). 

Stationary phase in the y'-variables gives the critical point equation 

\q{y)-q{y')\ \q{y)-q{y'\) 

The non-degeneracy follows from the identity 

\detWl,[y'9-i^{y,9)-\q{y')-q{y)\]\ = \detVl,My')-q{y)\)\ 



det Vy (Ty , -— — --t) 

\q{y')-q{y)\ 



>CaQ(z/j,,r(g(|/),g(y))) (1 + — — — - 



qiy)\ 

where Cgn > 0. The Hessian matrix is non-degenerate since the integral formula for c{fi) is 
microlocally cutoff away from grazing directions, ghost directions and unclean points and so, 
in particular, {i'y',r{q{y'),q{y))) > C{e) > 0. We denote the critical point by y'jjj,9) and 
let q'c = q{y'c{y,d))- Then, 

cf{jj,) = {2TTi2)-"+^[ f c^Mf-k(?/)-a^l+$(a:.afa).e)] ^^.o(g(^)^ gc^ IqJv) " gel) ^g^-^ 

JR"-! Jan \q{y) — q'^^ 

x(^,',r(g(y),g;))^ ■ (^1 + ^-^^^^^^^^^) ' Xii\y - y'W) daiy) d9 +i?(/x). 
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where, R{fi) = 0{jj, ") as /i — > oo. The critical set of the phase function \l/i(g(y), 6') = 

Crit(v^0:={(y,^)e9f]xM-^^ = (T,,,^|l^), ci,$ = 0, (T,, ^M^i|l) = d,,$}. 

The assumption |Sf(a)| = in Theorem [9] imphes that 

|Crit(^i)| =0. (95) 

The measure zero condition in ( l95l) follows by relating Crit(\l/i) to the Lagrange immersion of 
the critical set C$ = {t/^^ = 0} (resp. C^p^ = {d^^ipr = 0}), where friq, q', ^) '■= viji q, q', ^)- 
The Lagrangian immersions are given by 

i(j> ; C(j) — y 1 p 

{y, y', 9) ^ {y, dy$, y, -dy<^) = {y, 9; y, -dyip), de^ = 
and 

/ f \ , f f /rr (liy') - (liy) X /rj. q{y')-q{y)^\ , ^r/ n 

{y,y,u)^ ly,{Ty>,-—- -—),y,-{Ty,-—- -— , y = ttG {y,uj). 

V my) -q{y)\ my) -q{y)\ J 

It then follows that under the above identification, Crit(\l'i) corresponds to the set 

4w = {iy,v) e B*dn; (3{y,v) = f^Fiy,v)}- (96) 

By assumption, |SL'^J = and so |Crit(\l'i)| = also. Thus from flMl) and fl86l) it follows 
that cf (/i) = o(|U"""^) as /i — )■ oo and so, J^ e*^*/i^'^""-'^^cf'' (/i)d/i = o{t~'^) for |t| small. 

6.5. Multiple-link contribution to /''(t). The main result of this section is the compu- 
tation of the contribution of higher-order multilinks to /''(t). We prove here 

Lemma 6.5. Under the assumption that |Sp'(A)| = 0, 

I^.{t) = o{t~'') for all ] = 1, ..., M. 

Proof. We first compute the contribution to I^{t) given by double links with initial and 
terminal points pinned on the boundary. The computation for higher-order multilinks is 
carried out in the same way as for double links. In view of Lemma 16.21 the contribution of 
double links amounts to computing the small t asymptotics of 

/|W = E / xi'\s + t)([Ui':^]''{s + t)oF,{s)){q,q)dsdaiq). 

From fl82j) . wriing t =\z — g(|/')| the reflected parametrix 

\Uf}\\s + t, g(y), q{y')) = J^ t/g(. + t - r, q{y), z) o f/«(r, z, q{y')) dz 



= (27r/i)2" /^ /^„_, /^„_, /;° /;° ^^xny,y',r,.ur..,,.,r) B^'\q{y) , q{y'), z, r.co,, r,uj,) 
y<Xk{q{y') , r2UJ2)Xk{q{y) , riuiy'l'^r^'^ dridr2duiduj2 dz 
where, 

'if{y,y',riUJi,r2Uj2,z,T) := ri{q{y) - z,uji) + r2{z - q{y'),uj2) - (s + t - r)ri - rr2. 



(97) 
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In analogy with the single-hnk case, one apphes stationary phase in the spherical variables 
(wi,W2) e S^-^ X Sl-\ The result is that 

/f(t)= / (27r/.)2"e-^^*cf)(/i)ci/x, 



'0 
M(,,\.— II III I Jt^lri\q{y)-z\+r2\q{y')-z\-r-^{s-T)~r2T+s+<!>{q{y'),q{y),e)] A{k), 



Jn jR"-i J an Jan Jo Jo J i 
'xBfl{q{y), q{y'), z, rir{q{y), z),r2r{q{y'), z))x{ri-l) xi^\s+t)r'l-^r^-^dsdridr2da{y)da{y')dedz. 

(98) 

Now, stationary phase in (s,ri) gives the critical points ri = l,s = \q{y) — z\ + r = 
\q{y) — z\ + \q{y') — z*\ and so, the result is that 



c?H/i) = (27r/i)-W / / / e^^[l^(^)-^l+l'(^')-^l+*(^(^')'«(^)'^)U^5(g(|/),g(|/0,^,^) 
Jn iiR"-i Jan Jan 

X X^2\\qiy) -z\ + \qiy') - ^1) daiy)daiy')dedz, (99) 
with y4^o ^ C^{dQxdQxQxW^^^). We expand the z-dependent part of the phase function, 

^{z, y, y') = \q{y') -z\ + \z- q{y)\- z E n, (100) 

in normal coordinates {z',Zi) near the boundary, so that Zi = on dQ and Zi > in the 
interior of Q. Writing z = q{z') + ZiVz'i one makes a Taylor expansion around ^i = : 

^{z, y, y') = \q{y) - z\ + \q{y') - z\ 

= \q{y) - q{z')\ + \q{y') - q{z')\ + ((^.', r{q{y), q{z')) + r{q{y'), q{z')))) ■ z, + 0{zl). 

(101) 
Since the parametrix is constructed so that the operator wave front propagates along re- 
flected, non-tangential bicharacteristics, it follows that there exists a constant C{e) > Q such 
that 

(z/,,, r{q{y),q{z')) + r{q{y'),q{z'))) > Cie); {q{y) , q{y') , z) G supp4o(-,e). (102) 

By Taylor expansion of the amplitude, A^q, in ( l99l) around zi = and using the inequality 
(I102p . one integrates out the zi normal variable to the boundary in ( l99l) . This creates an 
extra power of /i~^ and the result is that 

4^)(/i) = {2nf^)-' /,^ /,^ 4„_, /,^ e^Mk(.)-.(^')l+k(.')-.(.')l+*(.(.'),.(.),^)l ai%{y), q{y'), q{z'), 6) 

y<X?My) - q{z')\ + \q{y') - q{z')\) da{y)da{y')d9dz' , 

(103) 
with a^o G C^{dVL x dVl x dVl x M"~^). One then applies stationary phase in the {y',z')- 
variables. Strict convexity ensures non-degeneracy of the Hessian. The critical point equa- 
tions are: 
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(T,,, r(g(y'), qim + dy,^q{y'), q{y), 9) = 0, 

(104) 
{T,:,r{q{z'),q{y'))) = -{T,:,r{q{z'),q{y))). 

The second set of equations in (11041) determines a locally unique critical point y'{y, 9) G dVt 
z'{y, 9) G dQ that is a reflection point for the double-link joining y G dQ and y'{y, 9) G dQ. 
The simplified formula is 

cf (/x) = (27r;u)-"-^ /" / e'^^'^y''^ 
Jan jk"-! 

X ciJ5(g(y), qiy'iy, 9)), z'{y, 9), 9) x^^\\q{y) - z'{y, 9)\ + \q{y'{y, 9)) - z'{y, 9)\) da{y) d9, 

(105) 

where, 

^2(2/, 9) = \qiy) - qiz'iy, 9))\ + \q{y'{y, 9)) - q{z\y, 9))\ + ^q{y'{y, 9)), q{y), 9). 

By another application of (186|) . one supplements the equations in (I104p with the critical 
point equations in {y, 9) : 

(T„ r(g(y), 2;')) + dy^{q{y'), q{y), 9) = 0, 

(106) 
. deHq{y),q{y'),9) = 

Let 

Cnt{^2) := {{y,0) e dQ xW^-';9 = (Ty,^^—^), d,$ = 0, (T,,, ^j^^) = ^^'*' 

|g(l/)-2;'| |g(y)-^| 

{n,,r{z\q{y'))) = -{T,,,r{z' ,q{y)))}. (107) 
Under the Lagrange immersion, t^^ : C<^^ — ^ L^, we have that 

V,(Crit(vI/2)) = sg;,), 

where 

41a) = {(l/>^) e fi*5fi;/32(y,r7) = «:p(y,r/)}, 

and so, the assumption that |S];,.%^J = implies and |Crit(\E'2)| = 0. Summation over the 
microlocal partition (ie. summing in k = 1, ...,M) and an application of fl86|) then implies 
that C2(yu) = o(/i-"-i) and so, ll{t) = J^ e^^*(27r/i)2"c2(/i)d/i = o(t-") as t ^ 0. 

The computation for arbitrary multlinks follows in the same way as for double links and 
this finishes the proof of Lemma 16.51 D 

Since l\t) = Y.f=i I-it) + R,{t) + R{t) where R{t) is locally smooth and i?^(t) = 0{e)t-'' 
near t = 0. Since e > in Lemma [6.11 can be taken arbitrarily small, by letting e — t- 0^ and 
applying the Fourier Tauberian theorem, this completes the proof of Theorem [9] under the 
assumption that dQ is strictly-convex. 
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6.6. WKB Case. In our application to QER, the semi-classical Fourier integral operators 
F{X) which arise are naturally represented in a somewhat different form with Schwartz 
kernels of WKB-type. We now modify the proof of Theorem [9] to apply more directly to 
operators of this form. The proof is very similar to the one given in Theorem [9] but simplifies 
slightly. 

Consider the special case where F{X) has Schwartz kernel 

F{\My),q{y')) = AVe^A^(''to).'^fe'))5(g(y),g(yO,A), 

where Biqiy),qiy'), \) r^x^oo EZoBMy),Qiy'))>^'' with B^ E C°°(9f] x dn) and supp 
Bk n Asuxdn = 0, k = 0,1, 2, .... The Fourier transform 

F{sMy),q{y')) = //x^e*'^(^(^(^)'^(^'))-) i?(g(y), g(y'),/i) ^Z^- 



To determine the leading (t + iO) ""-singularity of I^(t) = Tr f^U^(t + s)F{s)ds, one 
repeats the argument in Theorem [9l except that there is no 6'-integration at the end and one 
does not carry out stationary phase in y'. Rather, one directly computes that the single link 
contribution to the (t + iO) ""-singularity is given by f^ e*''*Ci(//)(i/i where, 

ci(/i) = (27r^)"- V/50 Lo e^A't^('^(^)-g(^'»-l'^(^)-^(^')ll ^»('^fa)'^fa2\ da{y)da{y') + 0(/i"-2) 

= Ian Ian Hf^Mv)^ ?(^')) [(Stt/x) Go{f^){q{y'), q{y))] da{y)da{y') + 0(/i"-2), 

(108) 
where, Go(/^) is the free Greens function in (I3T1) . 

The /c-link contributions for |fc| > 2 are computed in a similar way by inserting the 
reflection parametrix terms. The result is that 

cM= f f F{f,){q{y),q{y'))[{2iTf,)Go{fi)]'{q{y'),q{y))da{y)da{y') + 0{fi''-'). (109) 
Jan Jan 

As in the proof of Theorem[9l one shows that the critical point equations in {y, y') are equiv- 
alent to the fixed point equations f3^{y, rf) = Kpiy, rf); {y, rf) G B*dfi where {y, r]; Kpiy, v)) = 
{y,dy4'{q{y),q{y'))',y',—dy'ip{q{y),q{y')))- It follows from ( l86l) that in the sense of distribu- 
tions, 

I\t) = o(r") (110) 

under the measure zero assumption, |Sf(a)| = 0. 

6.7. Non-convex boundary. In the non-convex case, the canonical relation, Tj^an, of the 
boundary jumps operator is larger than graph of the billiard map, F/j. This is due to 
"spurious" geodesies that exit fi at a boundary point y G dQ, only to re-enter at another point 
y' G (9f2. Following [HZ], we say that {y,rj,y' ,rj') G B*dQ,xB*dQ with 77 = dy\q{y) — q{y')\ and 
T]' = —dy'\q{y) — q{y')\ is a spurious point of iyF((iV|^(A)) if the line q{y)q{y') leaves ^2. It is 
shown in [HZ] that these geodesies can be ignored for the purposes for spectral asymptotics 
of boundary traces of eigenfunctions. We briefly recall the argument here. Further details 
can be found in |HZ] section 11. 

Let 6 be a smooth, compactly-supported non-negative function on M" which vanishes on 
Q. Consider the metric 

gs = {1 + sb) gEuciidci^n; s e [0,6], 
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with S > small and the associated resolvent kernel Gs{z, z'; A) = (A^ — (A + iOy)~^{z, z') 
of the (yfs-Laplacian on M". The associated boundary jumps operator is 

[Nsf9^{X){q,q') = 2^G,(g,g', A); q,q' G dn., 



and one forms the averaged operator 

Nisw--= I x{s)[N,m\)ds, (111) 



where x^((0, 5)) with x ^ and /^ x = 1. The averaged operator, Nqq{\) satisfies two key 
properties (see [HZ] section 11): 

(OiV(A)< = ui, 

(112) 

The first condition follows from Greens formula since the metric gs stays Euclidean in Q 
and the second, by an integration by parts argument in s G [0, 6] (see jHZ] section 11). Now, 
let xT ^ C^{B*dQ) with xT = in an ^-neighbourhood of set of spurious directions and 
Xf = 1 outside a 2e-neighbourhood of the same set. Then, from (I112p 

= l^Ex,<x{NiX,)Opx,ixr)HX,)OpxMr)NiXjH^,ui^) + Oie) (113) 

= Nix) EA,<A(OpA,(xf )i^(A,)OpA,(xf )<., <) + Oie). 
In the second last line of (I113p . one uses boundary ergodicity and Cauchy-Schwartz to get 
the 0{e) error term. We replace F{Xj) by Opx.{xf')F{Xj)Opx.{xf') i^i the Weyl sum and 
repeat the argument in the strictly convex case. Finally, we recall that by removing the 
cutoff Xe supported in e-tubes around the grazing and singular set Q U B^dQ, we obtain an 
0{e)-eTTOT. As in [Zj IMP] we can let e — ?• to complete the proof of Theorem O □ 

Remark: We may also consider the singularity at any other value t = to of 

J]e^*^^(F(A,)4,,<.). 
j 

It gives the growth rate of 

5^ e^*«^^(F(A,)n^^,<,) 

i:A,<A 

As an example, when -F(A) = A^|q(A) and under the non-looping assumption, it follows by 
the argument in Theorem [9] that for to 7^ 0, |Si?(;^)| = and so. 



^ E ^''"''KWhon) = 0(1) as A ^ oo. 



A,<A 
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7. Local Weyl law on H: calculation of Woo 

We use the local Weyl law for Fourier integral operators to calculate the limit state Uoo{,o). 
It is provisionally defined on C°°{B*H) by 

Woo(a) = ^lim — — "^ {Opx^{a)^x,\H,Vx,\H)LHH)- (114) 

We now show that the limit exists and we calculate it. We use the notation, 

Xj<X 

for the limit measure in the boundary case. This section of course depends on the choice 
of H, but the cutoff estimates above show that the singular and tangential sets, hence the 
convexity oi H, is irrelevant to the result. 

Lemma 7.1. The limit state Woo is well defined and is given by 

Woo (a) = Woo yjH ) (^Ph- 

Proof. We prove the Lemma by reducing the calculation to the local Weyl law along the 
boundary in Lemma L2 of |HZj . 
By (TTBI) . we have 

uZia) = ^lim -i^ 5^(ArJ[,(A,)*OpA,(a)iVj[,(A,)4^,w^^)^2(^^^ (115) 

^ ' A,<A 

By flT7|l . the limit in flllSp is the sum of a pseudo-differential Weyl law term and a Fourier 
integral term. By Egorov's theorem (see Proposition 14. 2p . the limit of the pseudo-differential 
term is 

^ ' A,<A 

To complete the proof, we need to calculate the limit of the Fourier integral term. We claim 
that 

Aj<A 

By Theorem [9l it suffices to show that the set of points where the broken billiard map Ph 
equals a power of /3 has measure zero. However, this is special case of our non-commutativity 
assumption: If (5h{C) = f^^iO then /3^/3iy(C) = f3^^^{(). Also, if /3h(/3^C) = f^'^if^X), then 
both equal /3'^+^(C). Hence, iS^lSniC) = l^aniOl^^ '^^ ^ set of positive measure if /3//(C) = /3^(C) 
on a set of positive measure and if, for some ^ 7^ 0, /3//(/3^C) = P^W^O 011 a set of positive 
measure. We claim that the first condition implies the second: let the first set be called B. 
Then [i^B fl B must have positive measure for some I by Poincare's recurrence theorem or 
by ergodicity of /3. 

D 
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8. Quantum ergodic restriction: Proof of Theorem H] 

We now have the ingredients to prove the main result. As mentioned in section ^ we 
prove Theorem m first for A-pseudodifferentials Op\{a) G Opx{Sj {T*H x [0, A,^^])) with supp 
a nrg^( TUB^dVlUQ ) = 0; fc = 1, 2. In §8. 41 we show that for a full-density of egenfunctions, 
L^-mass is not conentrated on the set TQ^^B^dflUTUQ) (see |HZj for an alternative proof of 
this fact using the boundary trace of the Neumann heat kernel). Using this fact, in §8.41 we 
extend the proof of Theorem |4] to arbitrary symbols a G 5*^/ {T*H x [0, Aq^]). We therefore 
assume here that the symbol a G 5*^/ {T*Hx [0, Aq ^]) satisfies the following support condition 
dH]) with arbitrarily small but fixed e > 0. 

We now reduce the proof of Theorem H] to Theorem [101 

8.1. Reduction to Theorem IIOL To prove quantum ergodicity for the eigenfunction re- 
strictions M^. = (fix^lnj = 1,2,... we put F(A) = N^^{X)*Opx{a)N^^{X) and use ([I7D to 
compute the variance for the eigenfunction restrictions, u\_. 



7TT Yl K'^PA,(a)v^AjH,<^A,|H) -Woo(0-(F))p 



N{X) 



Xi<X 



^Ei([^(^^-)-^-(^(^)H'<)i' + ^(i) 



A,<A 



Aj<A 

^ Mat E I {m\,;a,syujU^^iF))]u',^,u\^) p+-^ Yl l (i^2(A,;a,e)<.,<,) 1^+0(1). 

^ ^ Aj<A ^ ^ Xj<X 

(116) 
where the limiting state cjoo(o"(-f')) was determined in [|71 

The second line in (I116p follows from the special case of the boundary Weyl law of [HZj 
that says that j^ ^a <a II'"Aj||l2 ~ 1. From the same asymptotic formula, it follows by the 
Cauchy-Schwartz inequality that the first term in the last line of (I116p is bounded by 

^ Y ([i^i(Ai; a, e) - u^aiF))]* o [F,{X,; a, e) - u^ia{F))]ux,,ux^) + 0(1). (117) 

We know by Proposition 14.21 that for e > sufficiently small, 

F,{X;a,e)eOpx{S"f{T*dn)). 

By Theorem [9] and boundary quantum ergodicity for A-pseudodifferentials |HZ] . it follows 
that for £ > small, 

^5^([Fi(A,;a,e)-u;o,(a(F(A)))]*o[Fi(A,;a,e)-a;oo(^(F))K,,<.) = 0A(l). (118) 

^ ^ Xi<X 
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As for the second term in (11161) . by making the further decomposition F2{\j; a, e) = F2i{Xj] a, e) + 
F22{\j]a,e) with Tp^^^^^xie) = graph /3|^; k = 1,2, it follows from the inequality \a + 6p < 
2(a2 + 62) that 

^\2\iF2{\,;a,E)u%,ui^)\' (119) 

^ MX) ^ l(^2i(A,;a,e)n^^,«^;r +^ J] \{F22i\r,a,e)u\^,u%)\'. 

The remaining step is the proof of Theorem [TOl 

8.2. Proof of Theorem not We now prove a special case of Theorem [TOl where Kp = /3^ 

(resp. f3l[), and where F{X) = F2l(A; a,e) (resp. -F22(A; a,e). We re-state it for the sake of 
clarity in the special case. 

Theorem 8.1. Let F2i{X; a, e) be a semi-classical Fourier integral operator associated to the 
canonical relation P\j. Assume that {Ph,P) almost never commute. Then 

iH^MA)Ei(^-(^.'«'^K'<)i' = o- 

Similarly for -F22(A; a, e) . 

The proof of the general case is essentially the same. As discussed in the introduction, 
the proof is based on time- averaging as in (pTj) . 

Proof. In analogy with the boundary case treated in |HZ] we define the M-th order time- 
averaged operator: 

M 

{F2,iX;a,e))M ■■= j^J2i^!SWTF2kiX;a,e)[Nl^iX)r; k = l,2. 

From the boundary jump equation Nqq{X)u\ = u\, we have for k = 1,2, 

j^ J2 \{F2k{Xr,a,e)u',^,u',^) \' = ^ ^ |( (F2,(A,; a,£))M<.,<.) T- 

By the Cauchy-Schwartz inequality (combined with the local Weyl law ]^) ^a <a II'^a 11^ 
~A^oo 1), we have 

^ Y, \{{F2k{Xf,a,e))Mu\,u\;)\^ < ^ $^ ( {F2k{Xf,a,e))\, {F2u{Xf,a,e))MuX,u%)+o{l). 

(120) 



A -)■ ^ A A -)■ ^ A 



Remark: This use of the Cauchy-Schwartz inequality may seem rather un-natural. We are 
trying to prove that {F2k{Xj;a,e))Mu\.,u\,) — )■ as A^ — )■ oo, i.e. that F2k{Xj;a,e))MU\. 
is asymptotically orthogonal to u\,. Use of the Cauchy-Schwartz inequality converts the 
orthogonality condition to the condition that the average over the spectrum of the norms 
\\{F2k{Xj;a,e))Mul\\ -^ is of order 0{jj). 
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It suffices to bound the averaged matrix-elements on the right side of f ll20p . Writing out 
these terms exphcitly, one gets 

N{X} Z1aj<a(-^2/c(Aj; a, e))lj {F2k{Xj; a, e))MU^.,u^^) 

_ 1 sr^M v^ 

APN{X) Z^m,n=0 Z^Aj<A 

— m2 2^m^n,m,n=o Same + j^p 2^^^„^o same 

(121) 
In (11211) . we have broken up the double sum into the diagonal and off-diagonal terms. 

Due to the time-averaging over m, n in the sums in (I12ip it is useful to first slightly refine 
the microlocal cutoff condition in dH]). We first prove (I128p and Proposition 18.21 for this 
somewhat more restrictive class of symbols and then apply a density argument using the 
pointwise Weyl law argument in Lemma 18.51 to deal with the general case. 

We assume that \m — n\ < 2M; m, n e Z and fix £ > 0. We say that a link is e-transversal 
at the boundary if it makes an angle of at least e to the tangent plane at each point of 
intersection with dfl. We define the extremal link-lengths: 



LH{M,e):= max sup {\'K/3jjT-''^H{y,r]) - 7rr~''^H{y,v)\ 

|m-n|<2M(j^^^)gQM 



+ Yl k/3%(2/,r/) - vr/3'=-%(y,r/)| + \7c(]H{y,v) - vl ), (122) 



k=l 

where, ^2*^ := {{y, t]); dist( {y, rj), ^2M+i^^2M+i^'F) — ^}- Thus, Lh{M, e) is the maximum 
length of if-broken and e-transversal geodesies with at most 2M-reflections at the boundary, 
dQ. Similarily, we define 

L{p,e):= inf ^|7r/3'=(y,r^)-7r/3'=-i(y,r^)|. (123) 

Thus, L{p, e) is the mininum length of e-transversal geodesies in Q with p-reflections at the 
boundary, dfl. We define the constant 

C{M,e) = 2M + l + max{\p\;L{p,e) <LH{M,e) + l}. (124) 

To estimate the time-averaged matrix elements in (I12ip . in the next two sections, we 
require that the symbol a G Sj {T*H x (0,Ao]) satisfy the somewhat stronger support 
assumption 

dist (supp a, T!^n{S*c{M,s) U ^*ciM,e) U r) ) > £ > 0. (125) 

We deal with the case of general symbols at the end in section 18.31 
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8.2.1. Analysis of the diagonal sum. The first important point is that the diagonal sum is 
bounded by 0{jj). Indeed, the operators 

(126) 
in the diagonal terms are A-pseudodifferentials; i.e. for each n = 1, ...M, 

A^^r.{\j;a,e) e OpxiS'f{T*dn)). 

By the A-Egorov theorem for boundary A-pseudodifferential operators (see |HZ] Lemma 1.4), 
it follows that with ( := {y,ri), 

An,n{^,-a,e) = Opx, {-J^ >< cj[FU\a,e)F2k{\a,e)WC)\+On,e{\-')L^^L^ as A, ^ oo. 

(127) 
Substitution of (11271) in (I12ip together with the symbol compuation for F2fc(A; a,e)]k = 1,2 
in Lemma 14.31 and an application of the boundary local Weyl-law for A-pseudodifferential 
operators |HZ] gives 

liniA^oo ]v(A) Ea,<a( ^n,n(Aj;a,£) u\^,u\^) = JB*an^^) X |c^(i^2fc(A;a,£))(/3"(C)p7"^(C)rf2/rf^ 

= !B*anT\0 X HF2k{ka.emC)\^dyd^ 

= lB*H(^i^^'^'^^)Pdnis,T)dsdT. 

(128) 
since /3 is symplectic. Summing over the M terms and dividing by jp gives the following 

formula for the diagonal sum 
1 ^^ / 1 \ I f 

n=l ~^°° y '^ '' Aj<A J Jb*H 

(129) 

8.2.2. Analysis of the off-diagonal sum. The next step is to apply Theorem [9] to the av- 
eraged operator Fm(A) in (1231) . The non A pseudodifferential pieces of this operator can 
be written as a sum of the A-FIO's Fm,n(A;a,e) where m ^ n and iyF((Fm,n(A; a, e)) C 
graph [(/3^)*(/^™')(/^'^)(/^^)] (^se (I130p ). One can locally write these canonical graphs in terms 
of generating functions and so the operators -Fm,n(A) can be written as a sum of A-FIO's with 
phase functions of the form fl5^ . Thus, Theorem [H] applies to these operators. 

We first observe that we can remove the outer factors of Nq^ since they have eigenvalue 
1 on M^.. For a G S^l (T*H x (0, Aq]) satisfying the support assumption in (11251) . we define 

F^,„(A,;a.£) := F2fc(A,; a,^)* [iV|^(A,)]™ [iV|^(A,)*]"F2fc(A,; a,^), m ^ n. (130) 

The remaining step in Theorem H] is the following 
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Proposition 8.2. Let a e Slf{T*H x (0, Aq]) satisfy the support assumption in U25\) . 
Then, under the quantitative almost commutativity assumption in Definition it follows 
that 

M 

Y, Y.^F^^^{\y,a,e)u%,u\) = Q. (131) 

■n,m=0;m^n Aj<A 



'vz^timw 



Proof. In the sum in f ll3ip . the quantum observables are non-pseudo-differential A-Fourier 
integral operators, i.e. their canonical relations are graphs of non-identity maps or corre- 
spondences. We apply the Fourier integral operator local Weyl law of Theorem [9] and the 
quantitative almost nowhere commuting condition in ([3]) to prove f ll3ip . 

For m 7^ ra, we denote the canonical relation of F™,„(Aj; a, e) by Tp^^^^^.-s)* [7v|g{A,)]™ [JV|^(a,)*]" F^kiXj-^s) ^ 
By the composition theorem for semi-classical Fourier integral operators, this canonical re- 
lation is the union of the two branches {k = 1,2), 

= graph [(/3^)*(/3™)(/3")*(/3^)] C B*dQ x B*dQ, 

where f3 : B*dVL — )■ B*dVL is the usual billiard map of the domain Vt and /3^ : C\£'^^(E) — )■ 
B*dQ are the branches of billiard maps defined in ( 1T3|) and in Lemma [4.11 Equivalently, the 
canonical relation is the graph of 

(/3|)"'/3"^-"/3^; A: = 1,2. 

Note that there are two factors of Ng^Ng^, which accounts for the two factors of {/3fj). 
By Theorem |9l the fixed point set whose measure we need to calculate is the fixed point 

set 



pm—n nk 



{(3'k)-'r-''l3'k{y,v) = iy',v'),G-^iy',v') = iy,v), r = V^^_„(g(y), g(y')), 



(132) 



where ipm-n is the phase function of Fm,n(A; a, e). Clearly, ipm-niQiy) , liv')) is the length of 
the broken trajectory which starts at q in the direction rj + ^yl — \rj\'^Vy., breaks on if, then 
proceeds on a broken geodesic until it has made m — n reflections at the boundary. Finally, 
for the last two links the trajectory breaks again at H. Equivalently, 



i(A;a,e) 



Upez'^^p.m-"' 



(133) 



ca 



^p,m-n = {{y,v) e B*dn : /3"-"/3^(y,r7) = /S'^ /^^V , v)} ■ 

As this shows, Si?^_^(A;a,e) depends only on m — n. 
We compute 



1 1 



lim supjy^^oo limA->.oo 'n{X}m2 

< limSUPjv/^oo -g Z2\m+n\=l M ^ 

< ClimsupM^^jgE 



l^n,m=0;m^n l^Xj<X\ Fm,n['^j] O',^) U^.,U^, ) 



2A/ 



lim 



|m-n|=l^^"^A-^oo 



l^X.KX-^'rnn 



2M ,. 

|m-n|=l ^I'^A^oo 



7V(A) 



7V(A) 

J2x,<x(^mn{>^j;a,e)ui ut) 



{\j;a,e)ui u^) 



(134) 

The final step in the proof of Proposition 18.21 reduces to estimationt of the Weyl sum in 
the last line of (I134p in terms of the quantitative ANC condition in Definition [31 First, 
we need to define the relevant /ip^^-measures (see Definition [3]). To do this, we need to 



QUANTUM ERGODIC RESTRICTION 49 

further A-microlocalize the Fmn(A; a, e)-operators away from large-time iterates of the grazing 
and singular sets. When a G Sj (T*H x (0, Aq]) satisfies (I125p . by semiclassical wavefront 
calculus, 

WF;^{F2{X;a,e)) C ^a {n\k\<c(M,e){iy^v) e B*dn- dist{(3''{y,v), S*dn) > e, dist(/3'=(y,r7), S) > e}) , 

(135) 
where, ^a : B*dil — )■ B*dil x B*dil with iAiy,v) = {yiViViV)- Consider the closed set 
%{M,e) U ^c{M,e) ^ ^ with \Qc(M,e) ^ ^c{M,e) U T| = and the disjoint open set 

U^' ■= nik\<c(M,e){iy,v) e B*dn;distiP''iy,r]),S*dn) > e, dist(/3^(y,77), S) > e}. (136) 

By the C°° Urysohn lemma there exists x*^ ^ C^{T*dQ) with the property that Xe^iv-, v) = 
1 when {y,7]) G Gc^M,e) ^ ^*c{M,e) U T and xf (?/,^) = for (y,r/) G f/f . To simplify the 
writing in the following, we abuse notation somewhat and write x*^ instead of OpA(Xe^)- 
Clearly, /3'=(f/f ) C f/*i where for \k\ < C{M,e), 

Ue!k--= n {{y.v)^B*dn-di^i{P\y,v),S*dn)>e, dist(7r/3^(y,r/),S) >£}. 

\i\<C{M,e)-\k\ 

Let x^ ^ C'o" {.dVt) be a smooth cutoff equal to 1 on Q^rj^^ ^\ U S^(jy^ e) '-' "^ and equal to 
on f/j^. Since W F'y^{[Nll^{\)f) C graph /3^ and C{M,e) > 2M + 1 > |m - n| + 1, it then 
follows by semiclassical wave front calculus that 

F„,„(A;a,£) = ^m,„(A;a,£:) + Oe,m-n{^'°^)L-'^L^, where, 

F^,.(A; a,e) := (1 - xf )i^2(A; a,e)*(l - xf ) • iV|^(A)(l - xiU-„,|) ■ iV|^(A)(l - X.V„|-i) 

■ ■ ■ (1 - xf2)A^ln (A)(l - xfi)i^2(A; a, £)(1 - xf ) + a,™-n,M(A-°°)i2^L2; |m - n| < 2M. 

(137) 

Simlilarily, we define the microlocalized operator with Schwartz kernel 

G^(g(?/),g(2/')):=(l-xf)Go(A)(l-xf)(g(l/),g(z/')), (138) 

where Gq{X) is the free Greens operator in fl?Il) . For each p G Z, the composite operator 
Fmn(A; a, £:)[AG'q(A)P is a A-Fourier integral operator of order zero in the sense of fl51|) . 
Consequently, for {q{y),q{y')) G (9^2 x Sfi, 

Pp,m-n(9(l/),g(2/')) :=limsup—3Y|F™„(A;a,£)(g(i/),g(yO)| x |[AGg(A)]P(g(i/'), 9(l/))l < oo. 

A-s>oo A 

8.2.3. r/ie Q AN C- condition. We make the 

Definition 8.3. Given the mapu : 9i7x9fi — Ag^x^Q — ^ -B*(9f] withw{y,y') = {y,dy\q{y) — 
q{y')\), we define the measures (see Definition [7? 



We can now define the QANC condition in terms of the measures /x: 



p,m—n' 
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Definition 13. We say that {(3h, P) quantitatively almost never commute if 



limsup — y^ ^limsup/ip,fc(CCp,A: n U^'^) = 0. 

Remark: For general piecewise smooth domains, the measures fip ^ are given by formulas 
that are quite complicated because of multiple-reflection solutions to the commutator equa- 
tions. However, when dQ is convex, the measures n't j^ simplify and satisfy the estimates in 
Lemma [8.4( b) below. 

Lemma 8.4. Fix e > and let a G Sj (T*H x [Aq, oo)) satisfying the support condition in 
(a) For m,n & Tj with 1 < \m — n\ < 2M, 



lim 

A— >oo 



7TT ^{Fm,n{^j'^(^M^^)'^X,^'^X,) 



iV(A) 



Xi<\ 



p;L{p,£)<LH{\'m—n\,e)+l 



(139) 



where, the measures fil^rn-n ^'^^ defined in Definition \8.3[ 

(b) When dQ is convex and A C {B*dQ fl U^) is measurable, for all k with \k\ < 2M, 

f -1 -1 ^ 

with Lp{y,r]) := |7r/3P(y,r/) - nl3P-\y,r])\ + {n P^-^y , r]) - rtl3P-^{y,7])\. 
Proof. Part (a) follows from the formulas in section ESI Indeed by (11091) . 



limsup;,^^ j^ J2x,<x Fm,n{\] a, e)u\ u\) 



lim sup;^_ 



IX) A^A,<A ^m,n[^j', O, e)Uy.U^) 



N{X) 



N{X) 

< Ep-MP,e)<L^i\m-nle)+i Ico,,m^r. M\ [^ G'o{X)]P{q{y) , q{y'))\ X \Fm,n{^;a,e){q{y'),q{y))\ dq{y)dq{y') 

~ l^p;L{p.e)<LH{\m-n\,e)+l f^p,m-n{^^P,m-n) ■ 

(140) 



As for part (b), when d^ is convex, we prove the estimate in (I139p by substituting 
the explicit WKB formulas for the operators N^^{X) ^ and Nl^{X) [HZ] suitably A- 
micro localized away from the generalized grazing and singular sets. Since A fl WF^^{{1 — 
x5)^Iq (^)(1 - Xe!k^i) = for all k = ±1, ..., ±2M, from ([60]), the WKB formula for the 
(1 — xfi)^'!^ (A)(l — Xe\~i) [HZj and a standard stationary phase argument, it follows that 

F™,„(A; 5)(g(y), g(y')) -A^oo (2vrA)"f'e^^<-"(^(^)'^(^'))[5i°l„,,(g(l/), g(y')) 

+ A-ii?«„,,(g(y),g(y')) + ...]. (141) 

Here, for m > n, ipm-niliv)^ liv')) is the length of the locally unique if-broken, e-transversal 
m — n link joining q{y) and q{y'). For m < n it is minus the length of the link. In the same 
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way, for all p E Z satisfying L{p,e) < Lh^M^e) + 1 (see f ll24p ) it follows by wavefront 
calculus and repeated applications of stationary phase that 



[(l-xf)AGo(A)(l-xf)F(g(z/),?(z/0) 



cf)\G,{X){l - xf^^"" 

^,^^ (27rA) Ve^A^^(^(^)'''(^'))[Aj;)(g(2/),g(2/0) + X~'A^^i{q{y),q{y')) + ••■]• (142) 



In (11421) . for p > the phase ipp{q, q') is the length of the locally unique p-link joining qiy) 
and q{y'). For p < it is negative of the length. In the following, it will be useful to define 
the sum of these phase functions 

^m-n,p(g(i/),g(2/') := ^,^_„(g(y),g(y')) + ^p(g(2/'),g(l/))- (143) 

It is clear from flT23|l and flT22]l that for \m - n\ < 2M, 

inf |vl>^_„,p(g(y),g(yO)| > L(p,£) - Lh(M,£). (144) 

We write i/k = nf3^{y,ri); k = 1, ...,m — n — 1 for the intermediate reflection points of the 
geodesic joining yo = y and y^_n = y' and put L{q{yp),q{yp_2)) = |g(z/p) - g(l/p-i)| + 
|g(^p_i) — q{yp-2)\ (see Definition [3]). Using the convexity of dQ, it follows that for p> 2, 



n-l 



n-1 
2 



lA'^^hnftA n(ii'\\\ < r\P "^ ( kfafc+i)-gfafc)l A 2 TTP-l ( 

\^P,e{q{y):q{y ))\ ^ llfe=l \\q(y^ + ,)-qiy^)\ + \q(y^)-q{y,^,)\J ^ 1 lfc=l {{u^^y^^ ,r(q{yk),q{yk~l))+r(q{yk)Myk + l))) ^ 

< L{q{yp),q{yp_2)y'^Yll=\{^yk^r{q{yk),q{yk+i))y^ x Yll2\{uy^,r{q{yk),q{yk-i)) + r{q{yk),q{yk+i))y^ 

= 2-(P-^)L{q{yp),q{yp_2)r^ Ur=\{'^y.,r{q{yk),q{yk+,))y'. 

(145) 
In the last line of (I145p we have used that {uy^^,r{q{yk),q{yk-i))) = (z/j^^,r(g(yfc), g(yfc+i)))- 

Similarily, we let y^f = -nPuijj, v), y\L-n\+i = t^PiI (^"""^Pui^J, v) and yf = nP^PHiy, v),k = 
1, ..., |m — n|. It follows that with a constant CH,n > 0, 

\Bl^L,My)^iiy'))\ < 2-ii'"-"i-^ic7H,n^(?(i/&-n|),?(<-ni-2))-"^- (146) 

The estimate in (I146P follows in the same way as in (11451) noting that A^|^(A)(g(y), q{y'))- 
kernel has the additional {i'y,r{q{y),q{y'))-teTin in the numerator as compared with the 
Greens function GQ{X){q{y), qiy')). This accounts for the absence of the additional {i'y^,r{q{yk), q{yk-i))- 
terms in the denominator of (I146p . It then follows from Theorem [9] (see also the argument 



52 JOHN A. TOTH AND STEVE ZELDITCH 

in subsection I6.6p . the time-cutoff Lemma [6. II and (11441) that 

hmSUp^^^ Jfr^T.X.,<xFm,n{\]a,e)u\ U\) = hmSUP;,^,^ MaT EA,<A(^'n,n(^i' «'^)«L<.) 



< hmSUp^.,^ ^Ep;L(p,.)<LH(|m-n|,.)+l I /c«i(*_„,,)[^<^o(^)]''(?(^)>?(l/')) 



xF^,„(A;a,e)(g(?/'),?(2/)) dq{y)dq{y')\ 

< hmSUp^^^ C T.p-Mp,e)<L„i\m-n\,e)+l I Icrit{^m-n,p) 

xB^^l^My'),q{y))da{y)da{y')\ 



i^'^^-nAliy)My'))A^^}(q{y),q(y')) 



< ClEp;L(p.e)<L«(|m-n|,e)+l/crit(*,„ „) lA^iiM, Qiv')) B^^lnA^iv') ' 9(y))\ My)da{y') . 

(147) 
Substitution of the estimates (I145p and (11461) in (11471) gives 



hm sup;^_ 



(A) Z^\,<X^rn,n{^j',0',£)Ux.,U\.) 



NiX) 



^J2p;L(p,e)<L„(\m-n\,e)+l'^ P "i" Jc«t(*p,„_„) -^(^(^|m-n|)' '?(l/|m-n|-2)) ^ 

xL{q{yp), q{yp-2))~^ Ullli^yk^ r{q{yk), q{yk+i))y^da{y)da{y'). 



:i48) 



We rewrite (I147P by making the change of variables w : {y,y') ^-)■ {y , TfTy^iqiy) , q{y'))) 
{y,ri). The result is that 



lim sup;^_ 



N{X) Z^Aj<A ^m,n{^j] 0,} ^J'^Xp'^Xj) 
<Ep-L(p,e)<LH(\m-nle)+lIcO,m „nt/M ^p(2/, ^)"'^^m-n(/5/^(l/, r/))""^ X ULl7~\P''iy, V)) dyd?]. 

(149) 
This finishes the proof of part (b). 

D 



In view of the quantitative almost commutativity assumption, substitution of the bound 
(11391) in the last line of (I134p concludes the proof of the Proposition 18.21 D 



8.3. Completion of the proof of Theorem |4t general symbols. Given a G Sj {T*Hx 
(0, Ao]) we let xf e C^{T*dn) as above and write a = [(rj^)*xf ]a+ [1 - {r9n)*xf)]a. From 
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Proposition 18.21 and the diagonal step, it follows that 

limsup;,^^ IvJaJ Ea,<a l(C'PA,(a)y5A,|//, Vx,\h) - Jb-h (^Pandsdr]'^ 

-L.i/(1 -Xf)(s,r)a(s,r)pff,rfsrfr|2 (150) 

+3 I Ib*h ^f (^' ^)«(^' ^) Pandsdr f 
+3 lim sup;,^^ j^ I ( OpA, ( [(ri^c)*xf ]a) <^A, k, <^A> ) p. 
From Lemma 18.51 it follows that 

limsup-^|(Op,^([(rf^)*xf]a)^A>,¥^A>)P = Om(^), (151) 

and it is clear that 

2 

Xf (s, T)a{s, t) p^^dsdr = 0m(£). (152) 

B*H 

Finally, since dist( supp ( a - [(rj^)*xf ]« ), 'rmi%iM,e) ^ ^c{M,e) U T) ) > e, it follows from 
flT28l) and Proposition E2] that 

limsupA^^ jy|^^;,^,<;, |(OpA,(a - [(rj^)*xf ]a)(/?A, |h, V'aJ//) - /b*h(1 " xf )(«, ^)a(s, r)p|^rfsrfrp 
< jg^Eili (limsupA^oo ]vlAyEA,<A(OPA,(a- [(rj^)*xf]a)n,n(Aj;a,e)M^^.,Myj 

= JI L*/f (1 - Xf )a Pi^rfsdr + 0,(1) < § ||a|U^ + 0,(1) 

(153) 
as M — )■ oo. By the quantitative almost commutation condition in Definition [131 for the last 
term on the RHS of fll53p . limsup£_j.Q+ 0,(1) = o(l) as M — )• oo. The constant C = C{Q) > 
in the last line of f ll53p is uniform in e and M. Letting e — ;■ 0+ kills the first two terms in 
f ll5ip and fll52p and then taking M — )■ oo kills the last term in f ll53p . This completes the 
proof of Theorem H] D. 

Remark: As a check on the midline of the stadium, we note that in that case |Sm„| = 1 
for all m,n since /3^ commutes with the transmission map (3q^ (which equals a). Thus, the 
sets in the union in f ll33p are empty except when m — n = p and that one has full measure. 
It is simple to check that the limit is non-zero. 



8.4. Extensions and generalizations. We now observe that the proof so of Theorems [T] 
and m extend to general symbols and to Dirichlet boundary conditions. 
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8.4.1. Eigenfunction mass along tangential and singular sets. The extension of Theorems 
[1] and m in section 18.31 (see fllSip ) to general symbols a G S'^''^(T*H x (0, A^j"^] relies on 
showing that for a full-density of eigenfunctions, mass does not concentrate along tangential 
or singular sets. 

The point of this subsection is to establish this fact (see Lemma 18. 5p by using a pointwise 
microlcal Weyl law argument to show that there is no such mass concentraion on any closed 
set of measure zero. This method has the advantage of working whether or not Q has a 
boundary. An alternative method is to use the potential layer Nq^(X) to reduce the problem 
to dQ and then apply the argument in |HZj (see Lemmas 7.1, 9.2 and Appendix 12) which 
uses the Karamata Tauberian argument for the restriction of the Neumann heat kernel to 
the boundary diagonal. 

In the following we let {x', x„) G M"^^ x M denote normal coordinates near the interior 
hypersurface H C Q defined hj Q 3 x = qH{x')+Xni'x' and write (^', ^„) for the corresponding 
fiber coordinates defined by T*Q 3 C^ = C, +^n^x'- In the following, we denote the restriction 
operator to H hj •y^ : f \-^ fl^ and continue to write ( : B*H ^ S^fQ for the map 

C:(x',r)^(a;',e' + yr^^^-')- 

Let Z G B* H he a closed subset with \Z\ = 0. Then, by the C°° Urysohn lemma, we can 
choose Xsz{x',^) smooth and positive homogeneous of degree zero with the property that 

ioT{x',o'^s*Hn, 

' Xe,z{^', = 1 when dist {{x', 0, C (2) ) < e, 

(154) 
. XeM^',O = when dist((x',0,C(2))>2£. 

We define 

XsM^^O = Xe,zix',0 ■ xM, 
where, x('") ^ C'^(]R) equal to 1 near u = and the corresponding restriction xi^z ^ 
Co~(T*i7) given by 

X^,zW,e)--=XeM^';^'^U = 0). (155) 

Consider the operator Op{xe,z) ■ C"^(fi) -^ C"^(fi) and let U{t) = e**^^^ denote the Neu- 
mann wave operator on Q. We wish to compute the small time asymptotics for ^ . e^^^^\\Opx. {x^,z)''^\- \\'l'^(h) 
and then apply the usual Fourier Tauberian theorem to get large A asymptotics for the A- 
microlocalized Weyl sum ^ J2x,<x \\C>Px,{Xe,z)uxXLHHy 
For |t| < dist (x', 917), we have that 

^,,^(t,x',x') := E,e''^'\[Op{xf,zK]i^')? 

(156) 

= bH0piXe,z)Uit)0p{Xe,z)*7*H]it,x',x'), x' G H. 

Since H is an interior hypersurface, for |t| < dist (x', 9^2), bicharacteristics starting from 
x' G H do not intersect the boundary, dQ. Therefore, by propagation of singularities, modulo 
a C°°-error, we can substitute the interior parametrix Uo{t, x, y) = (27r)~" J^„ e*[^'^'^''^^^*l^lla(x, y, ^)d^ 
for U(t,x,y) in (11561) . By the usual composition calculus for Fourier integral and pseudodif- 
ferential operators, it follows from ( I156P that 

[/,,^(t,x',x') = (27r)-" / e-^*l«l|x.,2(x',0l'rfe + ---, (157) 
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where the dots denote terms with lower-order singularities at t — )■ 0. By making a polar 
coordinates decomposition in the fiber variables in fll57p and using that Xe,z is homogeneous 
of degree zero, it follows that as t — )• 0, 

V,^z{t,x',x') = (27r)-" (/^.^^ |xe,^(x',u;)prfu;) (t + ^0)-" + o(t-") 

(158) 
= (2vr)-" (L.,^ \xf,zW, a' T\x\ O dC) {t + ^0)-- + o{t-). 

Finally, an application of the Fourier Tauberian theorem, integration over x' E H on both 
sides of (11581) and using that Op{xe,z) = Op\^{xe,z) (since Xe,z is homogeneous of order 
zero) proves the following lemma: 

Lemma 8.5. Let Z C B*H he closed with \Z\ = and for arbitrary e > let xfz ^ 
C^{T*H) he the symhol m / f753]) . Then, 

li^^^P WxT E WOp^AxlzKWl.^H) = i^Trr^'d \xU^',a\'l-\x,adx'dA = 0{e). 

A->oo iV(A) ^^ 3^1 \Jb*h J 

Thus, for a full-density of eigenfunctions, mass does not concentrate on measure-zero 
closed subsets Z C B*H. 

In subsection 18.31 (see (I15ip ). we apply Lemma [531 with Z = Tg^i^^jn,! e) ^ ^c(m e) ^ '^) 
to extend Theorem H] to arbitrary symbols a G S^'^ . 

8.4.2. The Dirichlet Case. The proof in the Dirichlet case is very similar to the Neumann 
case. There are a few minor changes: the relevant density is pH{s,r]) = (1 — l?]]?)^''^ for 
Dirichlet and this follows from the result in |HZj for boundary ergodicity. Also, the Nq^- 
operator gets replaced by the operator with kernel Go{q,qH',\) and this operator has the 
same microlocal properties as Nq^. Finally, of course the boundary traces of Dirichlet 
eigenfunctions are v\ : dyip\. 

8.5. Nodal intersections with interior curves: Proof of Corollary [81. 

Proof. Let fi C M^ be a piecewise analytic planar domain. The corollary follows from 
Theorem 6 in |TZj provided we establish the following "goodness " condition on C dVt: 



Wx\\l^{C) 



Here, we write u'^ = ^\\c and Uy^'" = dy^(p\\c and a > is an arbitrary positive constant. 
Assume the eigenfunctions (p\ are L^-normalized in Vt so that j^ \(p\f'dx = 1. Then, by 
Theorem m it follows that for an ergodic sequence of the m^'s, 

lkAlU2(C) ~ 1 

as A — > oo. By standard sup-estimates, we also have that 

snp\d,ipxix)\ = OiX'/^). 

Corollary [8] then follows from (I159p since clearly A^^^ ^ e""^ for any a > 0. 

D 
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9. Generic interior hypersurfaces H c VI: Proof of Proposition [6] 

In this section, we prove the genericity of the condition of almost never commutativity 
stated in Proposition [61 As mentioned in the introduction, the non-commuting condition 
should be generic in a much wider sense. We present only a simple result here for the sake 
of brevity; we hope to give a more general one elsewhere. 

There are at least two standard ways to approach generic domains. Following [CPS[ IPSlj 
IPS2[ [S]), we may consider the Frechet space C°°{S^, Q) of smooth embeddings f : S^ ^ fi" 
(the interior of Q); one could relax the smoothness to obtain a Banach space C^{H,VL) for 
some k > 2. A second approach (see [U llFOj ) is to fix a hypersurface H (here, a curve) and to 
parameterize nearby hypersurfaces by flowing out H under vector fields along H, extended 
smoothly to a neighborhood of H. In either definition, an infinitesimal variation of H is 
given by a smooth vector field X along H. We say that a property is generic if it holds for 
a residual subset of the relevant space, i.e. a subset containing a countable intersection of 
dense open sets. 

The main idea is to prove that the sets CO^^p of Definition [3] are smooth hypersurfaces 
away from the corners. To do this we use a well-known transversality theorem. We recall 
that if F : M — )■ A^ is a C^ map of Banach manifolds, then x G M is a regular point if 
DFx : TxM — )■ Tp(^x)N is surjective. Also, y G A^ is a regular value if every point of F~^{y) 
is regular. 

Transversality theorem jAR[ [Ql [U] Let ^ : Ti x B -^ E be a C'' map of Banach manifolds 
with E.l-i separable. If is a regular value of ^, and if ^h = ^{H, ■) is a Fredholm map of 
index < k then the set {H eH :0 is a regular value of^n} is residual in %. 

We let X = S"^ and Y = Vt. The Fredholm condition is then trivial. We denote an 
embedding by / : S*^ — )■ i7 and let H = f{S^). The properties of / which concern us are 
properties only of the image H. The principal property is that the measure of COp^k is 
positive. We introduce the partially defined symplectic correspondence, 

<D,,p : n X B*dn ^ B*dn, <^kAH, C) = /3-^/3^*/3^/3^(C)- 

Here, /3^ denotes the (partially defined) inverse to /3^, which makes sense since (3^ is sym- 
plectic on its domain of definition. More precisely, we restrict ( to the domain of /3~^/9^*/3-'73|^, 
which depends (slightly) on H. Our aim is to prove that, for generic (i.e. a residual set of) 
H, the set {( G B*d^ : 15~''' j3^ j3^ [5^{Q = C} is a hypersurface (i.e. curve) away from the 
corner set in B*dQ and hence of Minkowski content zero. It suffices to prove that is a 
regular value of the map ^j^p{H,() — C? i-6- that D{^jp{HX) — C) is surjective for each 
{H,C) such that <^j,p{H,C) -'C = 0. If %,p(if,C)C then b^<^j^p{H,C) ■ T^B*dQ -^ T^B*dQ. 
Surjectivity of 

/ - D(^^j^p{H, C) ■■ T^B*dn -^ T^B*dn 

holds if and only if 1 is not an eigenvalue of this kind of 'Poincare map' for generic H except 
for a possible curve in B*dQ. Note that if COpj had positive measure, then all eigenvalues of 
this Poincare map would equal one at any point of density {y, rf) of this set. In ^PSll IPS2| [S] 
the somewhat analogous result is proved that for generic domains, the spectrum of the 
Poincare map for every periodic reflecting ray omits 1 from its spectrum. It is possible that 
the multi-jet transversality approach of that paper could be adapted to our problem. Here 
we opt for a simpler proof when the billiard flow is hyperbolic. 
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Proof. We tacitly view the equation (3^/3^(11, v) ~ P'' PniyiV) = as written in local coordi- 
nates in B*dfi. We henceforth write (3'^ = (3h and suppress the fact that (3h is double- valued 
since we can separately consider the four branches of ^j,p- We would like to show that the 
derivative in H of the vector valued function PuP'^iy, 'n)~ P^ Puijli v) spans the tangent space 
TfSj(Sjj{y^ri)B*dQ at each point where the equation holds. 

The variation Snl^Hiv, v) is defined as follows: Let H^ be a smooth curve of hypersurfaces 
through Hq = H with variational vector field X along H. Then, 

It is a tangent vector to a curve through Ph{C) depending on the vector field X on if defining 
the variation, i.e. SPniC) • ThT-L — t- Tpjj{^c)B*dVL. To see the dependence on X we give an 
explicit formula for Ph similar to that in §3.11 



(5H{y. V) = Eit]jiy, -v),y, -v) - {Pniy, v) - 4(l/, -v))Ciiy, -V)- (160) 



Here, tniy, —rj) is the time at which the reflected ray hits dO. and (hiv^ v) is defined in ( !38l) . 
This vector is smooth in the unit normal u j , ■,. The variation thus consists of two types 
of terms: ones where we differentiate the hitting times and one where we vary the normal 
direction: 

5Cf,(y,-r/) = 25(C-.,.^(,,„,)).,.^(,_,, (161) 

It is clear that both terms depend on (and only on) X j , _ . (we recall that q^iy, —1]) is the 
point where the trajectory from {y,ri) intersects H (see Definition 13. 4p . Thus, 6f3Ho{y,—v) 
depends on (and only on) the value of X at the point q^iy, —rf). 

The variation of the equation has the form 6(3H{f3^{y, —v)) = Df3''{6(3H{y, —v))- As just 
noted, SPuiy, —i]) depends on (and only on) the value of X at q^iy, —v)- O^ ^h^ other hand, 
6f3H{f3^{y, —v)) depends only on the point qff{(3^{y, —i])) where the trajectory intersects H 
after p bounces. We now prove that these points must be different. 

Lemma 9.1. Let Q be a hyperbolic planar billiard table. Let H G Q be a convex curve. Then 
for almost all {y,ri) G B*dQ, the set of < 2 points where the line segment y/rj hits H is 
disjoint from the set where the line segment P^(3H{y,v) hits H. 

Proof. The statement is equivalent to saying that the 'probability' in dVt x H oi pairs {x, q) 
such that the trajectory defined by {x, — |f5^) has the property that the point q where it 
first hits H is the same as the point q' where it hits H after p bounces from dQ. We may 
think of the orbit as beginning at q rather than x with the same initial direction — |f5^ 
and proceeding along the billiard flow of Q until it has bounced off dQ p times. We then 
terminate the trajectory when it hits H for the first time after hitting at the pth bounce 
point Xp. li q' = q for a set of positive measure in x, then g is a self-conjugate point in 
the sense that there exists a normal Jacobi field along the broken geodesic which vanishes 
at t = (at the point q) and also at a later time when it returns to q. In fact, it is a self- 
conjugate point in the stronger sense that there exists a positive measure set of directions 
in S*Q for which the billiard trajectory returns to q (and additionally so that there exist 
exactly p bounces before the return to q). This is trivially impossible for hyperbolic billiards. 
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since the (broken) Jacobi fields defined by the bilhard trajectories have no zeros and hence 
there are no self-conjugate points (see [W|). D 

Although we did not state it, the argument is valid in all dimensions. The assumption of 
hyperbolicity is more than necessary; for the Proposition, we only use that the set (g, r]) of 
directions of tangential conjugate points is a set of content zero in B*dQ. 

Henceforth we assume y ^ y'. Since X ranges over all vector fields, we may assume that 
X is zero at one of these points and non-zero at the other. Since -D/3~^ is symplectic (hence 
an isomorphism), the following Lemma is sufficient to prove surjectivity of the differential 
away from curve of {y,!]) E B*dQ where /3H{y,v){y:V)- 

Lemma 9.2. //dimfi = 2, then for all ( e B*dVt, {5(3H{q,v)} spans Tr^i^q^,^)B*Vt as the 
variation runs over vector fields on H, except in the case where (3h{(1,v) = {QjV)> '^■^- where 
q + t({q,r]) intersects H orthogonally. 

Proof. We choose the variations so that the variation vector fields X vanish at QHiP^PHoiO) 
and span T, ry„)^- We want to prove that DHf3''6Hi3Hiy,v) spans Tpkp^(^y^^-^B*dVL. It is 
helpful to use the following notation: a tangent vector v to B*dQ is said to be vertical 
(f G VB*Q if it is tangent to the fibers of the natural projection to Q; it is horizontal 
{v G HB*Q) if the momentum variable is fixed while the position variable changes. 

We first observe that q{y) + t({y, —rf) is a fixed ray independent of H. We note that the 
reflection from H depends only on the tangent line to H at the point of impact. In varying 
iJ, we can vary (and only vary) the tangent line at the point of impact and the distance 
along the ray to the tangent line. These parameters correspond to two variations of H: in 
the first, we fix the point of impact and deform H so that the tangent line at the point of 
impact varies; while in the second, we hold the tangent line fixed and move the point of 
impact along the ray qiy) + tC(|/, — r/). 

The possible reflected rays thus have the form 

F(t, ve) = q{y) + hC{y, -V) - 2(t2 - t){ug, ({y, -r]))ug, 

where ^2 = t2{t, 9) is the time at which the reflected ray hits dVt and where vq = (cos 6, sin 6). 
The reflected direction is = ({y, —rf) — 2(z/0, ^(y, —ri))^^, so 

I^H{t,e){y,-Ti) = {YX{y,-v)-'^{'^e,C{y,-v))'^e)- 

We note that the direction is independent of t. So for fixed 6, 

Here, 

dY dt2 

Further, 

OPh dY d 

We now check that these two vectors span TiSjji^y^ri)B*dVL. We have, ^{{vg, O^e) = i^^d^ C)^o+ 
{u0,Q^uq. The two terms are orthogonal, so the vector cannot vanish unless both coeffi- 
cients vanish, but that is also impossible by orthogonality. This tangent vector is vertical 
while the t derivative is horizontal so to complete the proof it suffices to consider the case 
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where ^C ~ 2(z/0, Oug = 0. The sum vanishes at ^ = if and only C = ^) the normal to H 
at the point of impact. But in that case, /3h{0 = (■ 

D 

10. NOTATIONAL APPENDIX 

(1) To notationally disinguish points of dQ and points of H, points of dQ are usually 
denoted by y or by q, while points of H are denoted by qn- By a slight abuse of 
notation, we also write y ^-^ q{y) € dQ for the parameter of a local parameterization, 
and s ^-)■ g//(s) & H for a local parameterization of H; we sometimes use s to denote 
a point of H. The dual coordinates are r] G T*{dQ) and r G T*{H). 

(2) 7rJ(^) is the projection onto the tangent space at g of ^ G 5**^2; we use the same 
notation for dfl and for H. 

(3) r{qH, q) = \'^"Z'^\ '- the direction vector joining q G dQ and qn G -fT. 

(4) 7/f resp. 7an are the restriction operators to H, resp. dQ. 

(5) M^: boundary traces of eigenfunctions (see table). Traces on H: u^ = ipxln] u^^"" = 

(6) 7(2/,r/) = yr^H^(seeiD). 

(7) The length of successive double links is Lpd/,-?]) := \Trf3^{y,ri)—nf3^^^{y,ri)\ + \nf3^~^{y,ri) — 
npP-\y,v)\. 



(8) {y, rf) is the single-link billiard trajectory starting at y G dVt with direction (^(y, rf) 
7] + a/1 — Ivl'^J^y and endpoint on dVl. 

(9) rj^: the transfer map from B*dQ to i?*if (Definition [2]). 

(10) r|^: the transfer map from B*H to fi*9fi (Q. 

(11) C: The domain of rj^ (Definition [MD; 



1^12) Q: The grazing set. Generalized grazing sets: Q* := Ukezf3'^{S*dQ) and ^^ 



uik\<Nf3'is*dn). 



N 



(13) S: The singular set. Generalized singular sets: S* = Ufcgz/3 (-B^^) and S 

U|fc|<Ar/3 (-B^fi). 

(14) From Definition O Definition 13.51 and ( l40l) : The graph of the transfer map is denoted 

by r,H = T„a u r^ 



H;2. 



(15) The e-regular set: Ue := {{y,T]) G B*dn; dist(rJ^(S*(9fi uguT)>e}. 

(16) Nqq{X) is the quantum transfer operator of (TT5|) . 

(17) Nqq{X) is the boundary integral operator of (ITSl) . 

(18) /?// ( §4. H and flT3|) ): the once broken transmission map (or correspondence) on B*dQ — > 
B*dfl. It has two branches (3^, {k = 1, 2) 

(19) F(A) := iVjf,(A)*OpA(a)iVf^(A) (see dH; 
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(20) From Lemma Wl] The canonical relation of Nqq{\)* Op{a) Nq^{\) equals 

^an ^dii ^^^ Ph Ph 

(21) F^AX): terms in the decomposition of F{X) (see Proposition 1^ 
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